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ABSTRACT 

Using the temperature and ionization calculated in our previous paper, we model 
the spectral evolution of SN 1987A. We find that the temperature evolution is directly 



^ . reflected in the time evolution of the lines. In particular, the IR-catastrophe is seen in 

Q , the metal lines as a transition from thermal to non-thermal excitation, most clearly in 

!>■ ; the [O 1] A A 6300, 6364 lines. The good agreement with observations clearly confirms 

I the predicted optical to IR-transition. Because the line emissivity is independent of 

' temperature in the non-thermal phase, this phase has a strong potential for estimating 

' the total mass of the most abundant elements. The hydrogen lines arise as a result of 

■ recombinations following ionizations in the Balmer continuum during the first ~ 500 

CN ■ 

. days, and as a result of non-thermal ionizations later. 

, The distribution of the different zones, and therefore the gamma-ray deposition, 

i is determined from the line profiles of the most important lines, where possible. We 

"fH I find that hydrogen extends into the core to ^ 700 km s~^. The hydrogen envelope 

Qh' has a density profile close to /) oc 

from 2000 - 5000 km s'^. The total mass of 
O ' hydrogen-rich gas is ~ 7.7 Mq, of which 2.2 Mq is mixed within 2000 km s^^. The 

• helium mass derived from the line fluxes is sensitive to assumptions about the degree 

■ of redistribution in the line. The mass of the helium dominated zone is consistent with 
J> . ~ 1.9 Mq, with a further ~ 3.9 Mq of helium residing in the hydrogen component. 

^ , Most of the oxygen-rich gas is confined to 400 - 2000 km s~ , with a total mass of 

^1 ~ 1.9 Mq. Because of uncertainties in the modeling of the non-thermal excitation of 

■ " " ' the [O I] lines, the uncertainty in the oxygen mass is considerable. In addition, masses 

of nitrogen, neon, magnesium, iron and nickel are estimated. 

The dominant contribution to the line luminosity often originates in a different 
zone from where most of the newly synthesized material resides. This applies to e.g. 
carbon, calcium and iron. The [C I] lines, mainly arising in the helium zone, indicate 
a substantially lower abundance of carbon mixed with helium than stellar evolution 
models give, and a more extended zone with CNO processed gas is indicated. The 
[Fe II] lines have in most phases a strong contribution from primordial iron, and at 
t ^ 600 — 800 days this component dominates the [Fe II] lines. The wings of the [Fe 
II] lines may therefore come from primordial iron, rather than synthesized iron mixed 
to high velocity. Lines from ions with low ionization potential indicate that the UV 
field below at least 1600 A is severely quenched by dust absorption and resonance 
scattering. 
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Subject headings: abundances - line formation - nucleosynthesis - supernovae: general 
- supernovae: individual (SN 1987A) - stars: evolution - stars: interiors 

1. INTRODUCTION 

In Kozma & Pransson (1997; hereafter Paper I) we calculated the temperature and ionization 
in the ejecta of SN 1987A, which determine the line and continuum emission. Because one of 
the most important goals is to derive abundances and masses of the synthesized elements, these 
parameters are crucial for understanding the conditions in the ejecta. As will be shown in this 
paper, many lines have strong contributions from several abundance zones, and the emission may 
in fact be dominated by zones where only a small fraction of the mass of the element resides. In 
Paper I, it was shown that the temperature of the different abundance regions differ considerably, 
which puts an analysis based on a uniform temperature into question. In addition, the emission 
from a given zone is influenced strongly by the composition. Even a trace amount of an efficient 
cooler, like Ca II, can quench the emission from other lines. Finally, there is an interaction between 
the different abundance zones, mediated by the radiation. An obvious example is the emission 
from the hydrogen lines, which before ~ 400 days is powered mainly by UV emission from the 
other regions absorbed by the Balmer continuum. 

In spite of these complications, there has been important progress based on more limited 
forms of analysis. In a series of papers Li, McCray &; Xu, with coworkers, (Li & McCray 1992, 
1993, 1995, Li, McCray & Sunyaev 1993, Xu et al. 1992) analyze the most important emission 
lines from SN 1987A. Additional discussions are found in Kozma &; Pransson (1992: hereafter 
KP92), Wang et al. 1996, and Chugai et al. (1997). 

Prom the time evolution of the hydrogen lines Xu et al. (1992) and KP92 find that for times 
earlier than ~ 400 days photoionization from the n = 2 level dominates the ionization of hydrogen. 
In KP92 it was shown that the UV-emission from the ejecta, emitted as a result of the gamma-ray 
thermalization, could provide the necessary source for the Balmer continuum. Recently, Chugai 
et al. (1997) have modeled the time evolution of Ha, using a time-dependent model for the 
hydrogen-rich gas, similar to that in this paper and in Pransson, Houck &; Kozma (1996). They 
find that they indeed get a good fit to the observations of Ha up to the last HST-observations at 
2870 days after explosion. 

An important clue to the conditions in the ejecta was obtained from an analysis of the [O 

I] A A 6300, 6364 lines, using the fact that the lines went from optically thick to thin. From this 
Spyromilio &: Pinto (1991) and Li &; McCray (1992) estimate the density of the oxygen emitting 
gas and the filling factor and temperature of the oxygen component. In a further paper the [Ca 

II] A A 7291,7324 and Ca II A A 8600 lines are studied by Li & McCray (1993). Their main 
conclusion is that the calcium lines do not originate from the newly synthesized calcium, but 
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from primordial calcium within the hydrogen-rich regions. A similar conclusion was reached by 
Fransson & Chevalier (1989) in the context of Type lb supernovae. 

In modeling the infrared emission lines of iron, cobalt and nickel, Li, McCray & Sunyaev 
(1993) conclude that the iron-rich component must have a filling factor ^ 0.30 in the core. 
However, they do not take the contribution of other composition regions into account. A large 
filling factor for the iron is also indicated from dynamical arguments (Basko 1994, Herant & Benz 
1992). Finally, based on modeling of the He I A 1.0830 /im and He I A 2.058 /xm lines Li & 
McCray (1995) estimate the helium mass as ~ 3 Mq plus a similar amount mixed with hydrogen. 

An interesting result by Chugai et al. (1997) is that they find that the intensity of the Fe II 
lines can be explained only if trapping of the positrons from ^^Ti is efficient in the iron-rich parts 
of the ejecta. Because Coulomb scattering is not efficient enough, this puts interesting constraints 
on the strength of the magnetic field. At earlier epochs trapping should be more efficient. 

In this paper we exploit our results in Paper I for the calculation of the line emission from SN 
1987A, using a time-dependent formalism with non-thermal processes included. 

Our basic model, assumptions, and physical processes included, as well as the explosion 
models we use, have been discussed in Paper I. A brief summary of this is given in section ^ In 
section ^ we discuss the calculation and importance of the line profile for constraining the density 
distribution, as well as the influence of dust absorption. Section |^ contains a detailed discussion 
of our results, while in section ^ we discuss the limitations of our modeling, and the sensitivity of 
our results to various assumptions and simplifications. In section ^ we discuss our results, and in 
section |^ we briefly summarize our main conclusions. 

2. SUMMARY OF THE MODEL 

Because all details of the model are given in Paper I, we only summarize the main features 
here. All important elements are included with the ionization stages important in this analysis. H 
I, He I, O I, Ca II, and Fe I-IV are treated as multilevel atoms. The equations of the ionization 
balance and individual level populations are solved time-dependently to include freeze-out 
effects. Also the energy equation is solved time-dependently, allowing for adiabatic cooling by 
the expansion. Radioactive decay of ^^Co, ^^Co and ^^Ti are responsible for the energy input 
in forms of gamma-rays and positrons. The thermalization of these is calculated by solving the 
Spencer-Fano equation for the abundances and electron fraction in the zone. 

As input models we use the lOH model (Woosley & Weaver 1986, Woosley 1988) and the 
llEl model (Nomoto, & Hashimoto 1988, Shigeyama, Nomoto, & Hashimoto 1988, Shigeyama 
& Nomoto 1990, Hashimoto et al. 1989) for the abundances of the elements within the different 
burning zones. The ejecta is divided into a number of concentric zones of different compositions, 
masses and filling factors. Although simplified, this mimics the mixing of the different burning 
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zones in the ejecta. Outside of the core we attach a hydrogen envelope with a density profile fixed 
to give agreement with the line profile of the Hq line (see section |4.1.1[ ). 



3. LINE PROFILES 

The best way of constraining the emissivity variation, and therefore the density distribution, 
comes from modeling of the line profiles of the relevant lines. For a spherically symmetric 
distribution, it can be shown that for an optically thin line, the emissivity, j{V), is given by 

(Fransson & Chevalier 1989). j{V) is proportional to the gamma-ray deposition, and if one can 
consider the gamma-ray source as central this can be transformed into 

The first of these approximations is reasonable for many lines. In particular, lines dominated by 
recombination, following non-thermal ionization, or lines dominated by non-thermal excitation 
belong to this category. Examples of the former class are the H I and He I lines, while the [O I] 
A A 6300, 6364 lines at stages later than ~ 1000 days belong to the latter. If one line dominates 
the cooling of a given region this will also be a good approximation. One example of this is the 
[Fe II] A 25.99 /^m fine from the Fe - He zone at t ^ 800 days. 

The assumption of a central gamma-ray source is doubtful at velocities outside of which there 
is a substantial amount of radioactive material, V^\, in this case V^i ~ 2000 km s^^. A general 
gamma-ray deposition function, D^, is discussed in KF92, defined so that the mean intensity, J^, 
can be written as 

_ L.,Dy 
^^-le^^iJLe' 

where i?core is the core radius. For such a general deposition function the density distribution can 
be obtained from 

Approximations for are discussed in KF92 for the case of uniform distribution of the radioactive 
material, and outside of the core. Within the core one finds ~ constant, and 

p{V)^y^^ forV<Vm. (5) 

At V ^ ^Ni one has oc (as for a central source), and equation (|2|) is recovered. The 

important result is here the correspondence between the density distribution and the derivative of 
the line profile with respect to velocity. 
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Observations indicate that dust was present in the supernova ejecta at ~ 530 days (Lucy et 
al. 1989, 1991). The onset of the dust formation is hkely to have taken place akeady at ~ 350 
days and was completed at ~ 600 days (Meikle et al. 1993, Whitelock et al. 1989). From the fact 
that the dust seemed to be opaque both in the optical and IR (Spyromilio et al. 1990, Haas et al. 
1990), it has been suggested that the dust is clumpy, with very high optical depth in each clump. 
The latter is required to explain the presence of line asymmetries even at day 1806 (Bouchet et 
al. 1996). Lucy et al. (1991) estimate a covering factor f^av ~ 0.4 of the opaque dust. A similar 
covering factor, /cov ~ 0.55, is estimated by Wooden et al. (1993) from their dust temperatures 
and luminosities, based on a maximum dust velocity of V^ust ~ 1800 km s~^. 

The presence of dust will decrease the escaping fluxes in lines and continua by a factor which 
conveniently can be parameterized by a covering factor, /cov) applicable inside a velocity Vdust- In 
this paper we treat this as a pure absorption process. The thermalization of the absorbed radiation 
is discussed in more detail in connection with the photometry in in a subsequent paper. Here we 
only summarize points relevant for this paper. To model the absorption by the dust we assume a 
model where optically thick clumps of dust form at 350 - 600 days. In the dust formation phase 
it is probable that the covering factor depends on time, while later it is likely to be constant, 
with the optically thick dust clumps expanding with the rest of the ejecta. We therefore assume 
a linear increase in the covering factor from 350 to 600 days. At that time the covering factor is 
/cov = 0.40. Dust is likely to form only in the metal-rich parts of the ejecta (Kozasa et al. 1991). 
We therefore assume that only emission inside the velocity Vdust is affected by the dust absorption. 
Here we take Vdust = 2000 km s~^, similar to the core velocity. 

The dust may have important consequences for the UV intensity in the ejecta. Multiple 
resonance scattering increases the path length of the photons, increasing the effective probability 
for absorption of the UV photons by the dust in the core. This may be very efficient in decreasing 
the ionizing intensity in the ejecta, and introduces a major uncertainty in the calculations. The 
effects of a lower ionizing flux will be discussed in section |5.3| . Besides absorbing the radiation the 
dust may also cool the gas efficiently. This was discussed in Paper I. 

In our models we calculate the line profile directly from the source function, S{r), and optical 
depth in the line, T(r), as function of radius. The intensity at a given velocity is in the red wing 

/■-Rdust , , 

h{V) = 27r[ (l_/^„^)5(r)(l-e-^W) r dr 

JvRo/Vo 

+ / 5(r)(l -e-^(^)) r dr], (6) 

^dust 

where Vq is the velocity at the maximum radius of the ejecta, Rq, and -Rdust = 
[(Vdust^)^ + (^-^o/Vb)^]^''^- We have here neglected scattering of the background continuum, 
responsible for the P-Cygni absorption. More general cases, including this effect, are discussed in 
Pransson (1984). For an optically thin line without dust absorption equation (P) reduces to 

h{V) = 2ti / j{r) r dr, (7) 
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where j{r) is the emissivity. 

4. RESULTS 
4.1. Line Emission 

In order to compare our calculated line fluxes directly to observations we redden the calculated 
fluxes using Eb-v=0.06 from the Galaxy and Eb-v=O.W from the LMC (Sonneborn et al. 1996). 
The extinction curves are taken from Savage & Mathis (1979), and Fitzpatrick (1985). In addition, 
as explained in section ^ the internal dust absorption in the ejecta is taken into account for the 
emission from the core region. 

4.1.1. Hydrogen 

The hydrogen lines are discussed in detail in Xu et al. (1992) and KF92, and we will therefore 
be fairly brief in this respect. A major issue in connection with the hydrogen emission is the 
density distribution in the ejecta. As we discussed in section |^, the best way of constraining this is 
by the line profiles. Although a complete investigation of this question is a major, separate issue, 
we have made some experiments with different, simplified models. 

We divide the hydrogen distribution into a core contribution, most likely caused by mixing of 
hydrogen into the metal core, and an envelope component, which is probably relatively undisturbed 
by the mixing. Of these, the core component is especially uncertain. In the simulations by Fryxell, 
Miiller, & Arnett (1991) they find a hydrogen mass of ~ 1 M0, inside of 2000 km s-^ while 
Herant & Benz (1992) find ~ 2 Mq. Typical filhng factors are 20 — 40 %. We therefore treat the 
hydrogen core mass as a free parameter. For the envelope we use the density distribution of the 
Shigeyama & Nomoto (1990) 14E1 model, as well as our own, parameterized models. 

In Figure |l] we show the observed Ha line profile at 804 days, taken from Phillips et al. 
(1990), together with our models. Because only the envelope contributes at velocities greater than 
Vcore from the line center, we first discuss the effects of this on the line wings. The flux in the 
blue wing is complicated by the fact that Ha is optically thick. Photons from the background 
(i.e., continuum and other weak lines) are therefore scattered by Ha, and a P-Cygni absorption 
will result. Although we take optical depth effects into account for the line itself, we do not 
include scattering by other lines and continua. This most likely explains the asymmetry of the 
blue and red wings of Ha, and for this reason we concentrate on the red wing in this discussion. 
Dust scattering is, however, included in the emission from the core and behind it. The bump at 
933 km s~^ is caused by an unresolved [N II] A 6583 line from the circumstellar ring, as can be 
verified from high resolution observations. A weaker bump at 682 km s^^ corresponds to the [N 
II] A 6548 component of this doublet. 
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From Figure |T| it is clear that the 14E1 model gives too low a flux in the red wing. The reason 
is that the density profile in the envelope of this model is too steep. The llEl model, with a lower 
envelope mass, has an even steeper density gradient, and therefore does not improve the fit. For 
this reason we have simply parameterized the envelope distribution hy p = po (Vq/V)"', where po 
is the density at velocity Vq. po is determined by the total mass of the hydrogen envelope in the 
model. The exact comparison is somewhat sensitive to the continuum level assumed. As equation 
dl]) shows, the derivative of the line profile is, however, more relevant than the absolute level, and 
we find that a p oc V~'^ model, having a similar slope in the wing as the observation, is for this 
reason the most satisfactory model. The resulting line profile at ~ 800 days is shown in Figure 
H, and the fiatter density distribution improves the fit to the red wing considerably. However, 
between 5000 - 6000 km s"^ the p (X V model has a derivative that is too steep, and the 14E1 
model is closer to the observed slope. The mass within this velocity region is in the 14E1 model 
0.24 Mq, while in the p oc V~'^ model the mass between 5000 - 6000 km s"^ is 1.8 Mq. The 
maximum mass in this region should therefore be close to that of the 14E1 model, and ^ 0.5 M©. 
Our favored value of the total mass outside Vcore = 2000 km s~^ is therefore ~ 5.5 Mq. 

Having fixed the envelope component, we now add a core contribution, which gives the 
additional flux needed inside V^ore- The fact that the line profile is clearly peaked {dIu/d\V\ > 0) 
to ^ 700 km s~^ shows that hydrogen is present at least to this velocity. We find that a fairly 
uniform distribution of mass ~ 2 Mq in the core between 1000 — 2000 km s~^ gives an acceptable 
fit to the profile. Between 700 — 1000 km s~^ we need ~ 0.2 Mq. The mass within the core is 
therefore ~ 2.2 Mq, and the total mass of hydrogen-rich gas in the ejecta 5.5 + 2.2 = 7. 7 Mq. Of 
this ~ 3.9 M0 is hydrogen, while most of the rest is helium. 

With the hydrogen density distribution determined by the line profile, we discuss the light 
curves for the model with MH(core) ~ 2.2 Mq and penv oc V~'^. In Figures ^ and |3| we show the 
time evolution for Ha, Paa, Bra, Br7, the Ha/H/3 ratio, the Balmer, and Paschen continua, and 
the H 9^7, and H 7— >6 transitions. The solid line is the total calculated line fiux, using a full 
hydrogen atom with all nl-states up to n = 20 included (see appendix in Paper I) , which should be 
compared to the observations. In the figures we also show the contributions from the core regions 
(dotted line), and the hydrogen envelope (dashed line). Of these, the contribution from the core 
regions dominates up to ~ 700 — 900 days, after which the envelope contributes most. This is a 
result of freeze-out, which is most pronounced, in the outer, low density, envelope regions. 

In our calculations Ha is optically thick in the core up to ^ 1300 days. Ha becomes optically 
thin in the outer parts of the envelope already at ~ 550 days, and the entire envelope is thin at ~ 
850 days. H/3 is optically thick in the core regions up to day ~ 1000, while in the outer regions 
of the envelope H/3 becomes optically thin after ~ 350 days. The inner envelope regions become 
optically thin at ~ 700 days. The total Ha/H/3 ratio, as well as the ratios from the core and 
envelope regions, are shown in Figure |2[ When the Balmer lines are optically thick, i.e. Case 
C, H/3 splits into one Paa and one Ha photon. As H/3 becomes thin the Ha/H/3 ratio settles at 
a value of ~ 3.7 at 1200 days, slowly increasing to ~ 4.3 at 2000 days. Taking reddening into 
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account, the ratio should be multiphed by a factor of 1.18, resulting in an observed Ha/H/J ~ 4.3 
at 1200 days, as seen in Figure ^. 

Martin (1988) and Storey & Hummer (1995) discuss hydrogen recombination at low 
temperature. Martin does not take collisions into account, i.e. recombination in the low density 
limit. He compares calculations using an n- method (the different Z-states are populated according 
to statistical weights), and an n/-method where he treats the Z-substates explicitly. In Case B 
he finds that at 500 K the Ha/H/3-ratio is 2.5 and 4.2 for the n and n^-methods respectively, 
showing the importance of treating the Z-states individually. The assumption of complete Z-mixing 
(the n-method), however, becomes better for higher densities. Storey & Hummer include both 
individual /-states, as well as collisional processes. For an electron density of 10^ cm~^, a 
temperature of 500 K, and Case B they find an Ha/H/3-ratio of 3.68. The collisional processes 
makes the n-method a better approximation and reduces the ratio. We find in our models that 
/-mixing is accurate up to ~ 500 days. At later time the decreasing density causes increasing 
deviations from this approximation. At 1200 days this gives differences by ^ 50% for lines like 
Ha and Paa. At the same epoch the Ha/H/5 ratio is a factor of about two higher in the model 
with individual /-states. Therefore, at least at late time an accurate calculation requires inclusion 
of these effects. In our calculations at 2000 days Case B is valid, but our temperature, as well 
as the electron density, is lower than in the example from Storey & Hummer, resulting in our 
somewhat higher ratio. The effect of decreasing temperature and electron density is to increase 
the Ha/H/3-ratio. 

Photoionization from n = 1 is never a dominant process. Instead, for t ^ 500 days 
photoionization in the Balmer continuum is the most important source of ionization, as was found 
already in KF92 and Xu et al. (1992). Here we demonstrate this in a more realistic context, with 
a more self-consistent treatment of the UV-field, temperature and ionization. At t ^ 500 days 
non-thermal ionization from the ground state dominates photoionization from excited levels, in 
agreement with KF92. 

For hydrogen we include non-thermal excitations up to n = 4. For n = 2 the contribution from 
non-thermal excitations is larger than direct recombination, except for the outermost envelope 
regions at later times. For n = 3 and n = 4 recombination dominates during most of the evolution. 
Later than ~ 1000 days, however, for n = 3 the two contributions are of the same order, while 
for the n = 4 level recombinations always dominate. Within the core two-photon emission always 
dominates the de-population of n = 2 over Lya emission. In the envelope, however, Lya emission 
exceeds two-photon emission. 



4.1.2. Helium 

The regions contributing to the He I lines (see Fig. ^) are naturally the helium, but also 
the iron and hydrogen regions. In our calculations we find that the fraction of the deposited 
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non-thermal energy going into these three regions are ~ 8%, ~ 5%, and ~ 60%, respectively, at 
around 500 days, in accordance with Table 3 in Paper I when allowance for the positron input to 
the Fe - He zone is made. For determination of the helium mass the distribution in velocity is 
extremely important, since this directly determines the gamma-ray deposition. 



K., AM(V) 



A relatively large mass, AM(y), at high velocity can therefore be difficult to detect. This is 
further discussed in Houck &; Fransson (1996) for SN 1993J. 

In the same way as for Ha, we constrain the helium distribution from the line profiles. Of 
the He I lines the A 2.058 /im line is best suited for this type of analysis, because of the relative 
absence of blending with other lines (Fig. ^). For this purpose we take the observations by Meikle 
et al. (1993) at 574 and 695 days, when the line is optically thin (see below). Although the line is 
relatively free from blends, there may still be weak lines superimposed. Also, the signal to noise is 
limited. The helium density distribution based on this line fit is therefore uncertain, especially for 
V ^ 3000 km s~^. Because of the change in flux during this time interval, we normalize the line 
profiles to the peak flux of the line. The total line fluxes agree well with the observations, as can be 
seen from Figure |^. We flnd that the helium density is relatively flat between 1500 — 3500 km s^^, 
and then falling above ~ 4000 km s~^. Therefore, most of the mass is at velocities ^ 3000 km s~^. 
The continued rise of the line proflle even inside of 1000 km shows that there is a substantial 
amount of helium close to the center. The extension to ^ 4000 km s~^ is in contrast to the model 
by Li & McCray (1995), who use Fmax = 2500 km s~^. Because most of the helium is outside of 
the core, only the fraction in the Fe - He zone and the helium core fraction are affected by dust 
absorption. In our calculations we have 0.6 Mq of He within 2000 km s~^ and 1.4 Mq between 
2000 and 4000 km s~^, i.e. a total helium zone mass of 2.0 Mq. This is the total mass in the 
helium zone, of which ~ 1.9 Mq is pure helium. In addition to this, we have ~ 3.9 Mq of helium 
from the hydrogen-rich regions. 

Figure]^ shows the A 1.083 fim and A 2.058 fim light curves. The He I A 2.058 fim line is 
dominated by emission from the helium regions at all times, and therefore the differences between 
the models are small. For the A 1.083 iim feature, on the other hand, we flnd that the contribution 
from [S I] A 1.0820 /im is important, and actually dominates the emission for model lOH, up to 
~ 700 days. At later epochs the He I line from the helium zone dominates the light curve. Later 
than ~ 1200 days the contributions from the hydrogen envelope takes over. From the hydrogen 
envelope, most of the emission is due to He I, but the contribution from Pa7 increases with time. 

The flux of He I A 2.058 /xm is sensitive to the treatment of the continuum destruction 
probability of the A 584 line. The large abundance of carbon in the He - C zone can cause 
the A 584 line to be absorbed by C I. This process, and the competing processes of escape and 
branching to the A 2.058 yum line, has been discussed in KF92 and by Li &: McCray (1995). 
As was mentioned in the discussion of the continuum destruction in Paper I, we flnd that with 
the destruction probability from equation (32), in Paper I photoabsorption of the A 584 line is 
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important only for t ^ 300 days. The continuum destruction probability was in KF92 estimated to 
~7xl0~'' (r/5000K)^/2 [X(C)/0.018]. The reason why continuum absorption is at all important, 
given the branching probability of 1.1 x 10~^ between the A 2.058 /xm transition and the A 584 
transition, is that the optical depth of the A 2.058 ^m line is high. At 300 days the optical depth 
is in our standard model ~ 110, giving roughly equal probability of branching and absorption. 
After this epoch the optical depth rapidly decreases. Li & McCray (1995), on the other hand, find 
that up to ~ 500 days most A 584 photons are destroyed by this process in their He - C zone. 

To study the effect of the form of the continuum destruction probability, and therefore also 
partial versus complete redistribution, we show in Figure ^ the light curves for both the case of 
continuum absorption calculated assuming only absorption within the Doppler core (eq. [32], 
Paper I), and for the case of absorption dominated by the damping wings of the line (eq. [33], 
Paper I). From this figure we see that the different assumptions give total fluxes different by a 
factor of ~ 2.5 at 200 days, and a factor of ~ 5 for the emission from the He - C zone alone. The 
reason for the different factors is that continuum destruction is only important in the He - C 
zone, and not in the hydrogen zones. The factor of ~ 5 can be traced directly to the continuum 
destruction probability. For the damping parameter of the A 584 line, a ~ 2 x 10"'^, and a 
typical continuum to line opacity of A;c/A;l ~ 2 x 10~^, the Doppler case has a factor of ~ 8 
lower destruction probability compared to the damping case. In the damping case continuum 
destruction is consequently important up to 500 days, in accordance with Li & McCray, although 
the optical depth in the A 2.058 iim line is a factor 2-3 smaller. As we remarked in the discussion 
of the continuum destruction in Paper I, the more realistic case of partial redistribution gives a 
lower importance to the line wings, and a destruction probability closer to the Doppler case, and 
we therefore believe that this case gives the best approximation to the line flux. Chugai's (1987) 
partial redistribution approximation (eq. [34], Paper I) gives for fec/^L ~ 2 x 10~^ only a factor 
~ 1.4 higher destruction probability than the Doppler case. 

The effect of the continuum absorption can also be seen in a model, further discussed in 
next section, where we have replaced carbon in the helium zone by nitrogen, with an abundance 
of ~ 3.5 X 10^^. This model has a negligible continuum destruction of the A 584 photons, and 
therefore higher A 2.058 /im flux at early time. At 200 days the flux is a factor of 1.6 higher 
than in Figure ^ decreasing to 1.25 at 400 days, and 1.08 at 600 days. This again confirms that 
mass estimates based on the A 2.058 fim line are most reliable at ^ 600 days. As expected, 
the A 1.0830 /im line is not affected by this uncertainty. Instead, it is more sensitive to the 
temperature as well as the optical depth, and in addition blending with other lines. 

Observations by McGregor (1988), and Meikle et al. (1993) indicate a large optical depth 
in the two helium lines for the first couple of years, based on the asymmetry of the line profile, 
or rather the blue-shifted absorption troughs. In the spectra by McGregor the A 2.058 fim line 
is clearly asymmetric at 437 days, indicating an optical depth substantially larger than one. At 
574 days the observation by Meikle et al. is consistent with the line being either optically thick or 
thin, while at 695 days the trough has clearly disappeared. In our calculations we find an optical 
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depth larger than one in the helium zones up to ~ 700 days. The optical depths in the hydrogen 
core, and envelope are somewhat smaller, and the A 2.058 i-im line becomes optically thin in these 
regions at 650 and ~ 400 days respectively. 

For He I A 1.0830 ^m it is difficult to extract any information on the optical depth from 
the observations, due to the blending with other lines. In our calculations we find an optical 
depth greater than one in the helium region even for t ^ 2000 days. In the hydrogen regions the 
A 1.0830 fim line becomes optically thin at ~ 900 days. 

Li &: McCray (1995) find that in order to fit the A 1.0830 fim and A 2.058 fim emission they 
need ~ 3 Mq of nearly pure helium and ~ 11 of hydrogen mixed with primordial helium. 
They assume a filling factor of 0.30, but point out that their calculations are not sensitive to the 
choice of this parameter. Our helium mass is lower than Li & McCray 's. As we have discussed 
above, Li & McCray have a substantially larger destruction of the A 584 photons, and therefore a 
lower flux in the A 2.058 /um line than in our model. Consequently, at early time the contribution 
from their He - C zone is very low, and most of their flux originates at these epochs from the He 
— N zone. This probably explains their higher helium mass. 

In most other respects, however, our calculations agree well with those of Li &: McCray. 
In particular, we find the same distribution between the various contributions to the excitation 
of the 2p^P and 2p'^P levels. Non-thermal, direct excitation and recombination, following the 
non-thermal ionization, give in our models roughly equal contributions to the 2p^P level, each 
~ 40%. Earlier than ~ 600 days thermal excitations from the 2s^S gives an additional ~ 10 — 20%. 
The 2p^P level has a large contribution from recombinations, ~ 60% before 800 days, increasing to 
~ 97% at 1200 days. Because of the meta-stability of the 2s^S level, thermal collisions contribute 
~ 50% earlier than 800 days. At later epochs this contribution falls rapidly because of the 
adiabatic decrease of the temperature in the helium zone. 

Li & McCray find that in the hydrogen envelope the He I 2s ^S* state is depopulated by 
Penning ionizations. This is confirmed by our calculations. Photoionization from the ground state 
of He I is always unimportant. However, photoionization of excited levels is important for ionizing 
He I. Up to ~ 700 days the photoionization rate from excited levels in He I is somewhat higher 
than, or of the same order as, the non-thermal ionization rate. After ~ 700 days the non-thermal 
rate slowly becomes more important. The most important photoionization source is emission from 
lines in the UV, and especially the He I two-photon continua. 

4-. 1.3. Carbon and Nitrogen 

b igure I shows the [CI] A A 9824, 9850 light curve. Although the shape of our light curve 
is in agreement with observations, the model over-produces the line by a factor of ~ 10. Most 
of the contribution to the A A 9824, 9850 lines comes from the helium component, and the 
O - C region. The mass of the latter varies substantially between the llEl and lOH models. 
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0.10 Mq and 0.60 Mq, respectively. A possible ingredient in reducing the line strengths is the 
influence of CO. Liu &; Dalgarno (1995) find that while only a small fraction of the carbon goes 
into CO, the cooling of the gas is increased by up to an order of magnitude. The temperature 
consequently decreases at 500 days from ~ 3000 K without CO, to only ~ 1200 K including CO. 
This can easily decrease the [C I] emission by an order of magnitude. However, although the [C 
I] emission from the O - C region can be killed in this way, the emission from the helium region 
is more difficult to quench. As Liu & Dalgarno point out, CO is efficiently destroyed by He II, 
and little CO is expected to form in this region. One possibility is that we have under-estimated 
the photoionization flux above 11.26 eV in the model, which would explain the discrepancy. Our 
neglect of UV scattering argues against this. 

Because we have fixed the total gamma-ray deposition in the helium zone from the line profile 
and flux of He I A 2.058 ^m, the [C I] flux from this region should be fairly reliable. We therefore 
conclude that the most likely solution to the over-production of the [C I] line is that the carbon 
mass mixed with helium is lower than in the llEl and lOH models (^(C) ~ (1 — 2) x 10^^ in the 
helium region in both models). 

The enrichment of carbon in the helium shell occurs as a result of convection during the final 
helium shell burning phase (e.g., Arnett 1996). Both the time scale and the efficiency are, however, 
uncertain due to our limited understanding of the convection process. The amount of processed 
carbon mixed into the helium shell is consequently uncertain. To satisfy the observations, a 
decrease of the carbon mass in the helium region by a factor of 5 - 10 is required. 

To check the effect of more limited mixing of carbon into the helium shell, we have replaced 
the He - C zone by a zone with only hydrogen burning products, as given by the He - N zone in the 
lOH model. The most important difference is the high abundance of nitrogen, X(N) ~ 3.5 x 10~^, 
and low carbon and oxygen abundances, X(C) ^ 9.1 x 10"^ and X(0) ^ 5.9 x IQ-^, all by 
number. 

This model (Fig. ^) basically extinguishes the [CI] A A 9824, 9850 emission from the helium 
zone, as expected. The total flux is now close to that observed, especially taking the likely effects 
of the CO-cooling in the O - C zone into account. A test of this model is to check if the emission 
in lines of N I and N II is now compatible with the observations. The strongest of the nitrogen 
lines is the [N I] A A 10398, 10408 multiplet. Although not discussed previously, the CTIO spectra 
by Phillips et al. (1990) show a clear line at this wavelength in all spectra covering this wavelength 
region. On day 786 the model gives a reddening adjusted flux of ~ 6.5 x 10"^'^ erg cm~^ s~^. 
Including internal dust absorption this gives a flux of 3.9 x 10^^^ erg cm~^ s^^. On the same 
day the CTIO spectrum by Phillips et al. gives a flux of ~ 6.9 x 10~^^ erg cm~^ s~^ for the 1.04 
lim line, entirely consistent within the uncertainties of the model. The model strengths of [N II] 
A A 6548, 6583 are only ~ 1 — 2% of the Ha line. The blending of these with Ha will therefore 
effectively hide these lines. We therefore conclude that there is evidence from the observations for 
a more extended He - N zone, at the expense of the He - C zone. 
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Phillips k Williams (1991) argue that that the [CI] A A 8727/9824 + 9850 ratio is ^ 0.6 on 
day 589. On the same day we find that this ratio in the He - C zone is ~ 0.8, while in the O - C 
zone it is only ~ 0.3. With the 'standard' model the total ratio, which is essentially the ratio in 
the He - C zone, is in conflict with the observations. With the He - C zone replaced by a He - N 
zone, the total ratio is close to that in the O ~ C zone, and consistent with the observations. This 
provides some indirect support for our conclusions above. 



4-1 -4- Oxygen 

The [01] A A 6300, 6364 lines are of special importance for the analysis, because they are 
in an easily accessible, non-blended part of the spectrum. Therefore, a complete and accurate 
data set exists for these lines. In addition, oxygen is the most abundant of the metals, and a good 
probe of the progenitor evolution and its mass (e.g., Thielemann, Nomoto, &i Hashimoto 1996). 
For future reference we show in Figure ^ the levels and transitions included in our calculation. 

Starting with the lOH model, we show in Figure ^ the light curve of [01] A A 6300, 6364, 
together with the individual contributions from the different zones. Observations are taken from 
Danziger et al. (1991). The light curve can be divided into one epoch when thermal excitations 
of the level dominate, which lasts up to ~ 800 days, and one later epoch when non-thermal 
excitations dominate. This transition is, as explained in e.g. Fransson, Houck &; Kozma (1996) and 
as can be seen from a comparison with Figure 2 Paper I, intimately coupled to the temperature 
evolution of the core. As the oxygen-rich regions undergo a thermal instability and cool to ^ 300 
K, thermal excitation of the level effectively stops. The contributions from the hydrogen and 



helium-rich regions are, however, as shown in Figure |10|, significant in the thermal phase. In fact, 
up to ~ 30% of the [01] A A 6300, 6364 emission at 600 - 900 days comes from these components. 
Because of the adiabatic expansion and therefore falling temperatures of the hydrogen and 
helium-rich gas, the emission from both these decrease rapidly after ~ 800 days. This means that 
even if the non-thermal contribution from the core is small, it dominates after ~ 900 days. 

In terms of the qualitative evolution we find good agreement between model calculations and 
observations throughout the whole period, showing without doubt that the IR-catastrophe really 
has taken place in the oxygen-rich gas. Quantitatively, there are, however, some disagreements. 
While most of the thermal phase is well reproduced, at ~ 750 days the lOH model over-produces 
the [O I] luminosity by a factor of ~ 2. The llEl model shows good agreement with the whole 
thermal part of the light curve. 

The most serious flaw is in the non-thermal part. Although agreeing qualitatively with the 
observations, the level of the flat, non-thermal part of the curve is down by factors of 4 - 6 in the 
two models. As a consequence of this under-production, the transition from the thermal to the 
flat non-thermal light curve occurs 100 — 200 days later than is observed. In this context we note 
that the observed break coincides well with the calculated transition from the O - Si - S zone, 
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arguing for a larger contribution at late times from this. 

As an illustration of the sensitivity to the assumed mass of the oxygen region, we show in 
Figure the effect of varying the oxygen region mass by a factor of two from that of the llEl 
model, covering the range of M{0) = 0.95 — 3.80 Mq. Other parameters of the llEl model are 
kept constant. Here we first note that the non-thermal part is under-produced even in the highest 
mass model, so increasing the oxygen mass does not solve that problem. The thermal part of the 
light curve before ~ 500 days is best reproduced by the standard llEl model, and the two other 
extremes probably bracket the likely range of the oxygen zone mass. Between 500 - 900 days 
the difference between the models is small because of the large contribution from the helium and 
hydrogen regions. 

The total oxygen mass in the llEl and lOH models is roughly the same. In spite of this, 
the flux in the llEl model is only ~ 50% of that in the lOH model in the non-thermal part. In 
addition, although the mass of the O - C zone is only ~ 0.6 Mq, compared to ~ 1.2 Mq in the O 
- Si - S zone, the former contributes twice as much flux as the latter in the lOH model. 

To understand this, and the general level of the light curve we have to discuss an important 
technical point in the excitation to the level. During the non-thermal phase, in addition to 
direct excitation, this level receives an important contribution from the excited Ss'^S level, via 
the A 1641 line (Fig. P). Because excitation of the Ss'^S level is also non-thermal, mainly by 
recombination following non-thermal ionizations, the basic scenario is not changed, but it can 
make an important quantitative difference. Normally, the de-excitation of the 3s^S level is to the 
ground state, by the A 1302.2 — 1306.0 resonance multiplet, and the probability of a transition to 
the level in an individual transition is only A3s-^^_id^/ J2i=o ^3s-^^_3p. = 3.0 x 10~^. The optical 
depth of the resonance lines are, however, very large, increasing the effective lifetime of the 
level. The actual probability of a transition in the A 1641 line to ^D2, compared to is 



Pap 247r5i 

where the indices 1 and 2 refers to the 2p^ ^P ground state and the 35^5" state, respectively. 
Inserting atomic parameters, and using the O I density derived by Li &: McCray (1992), one finds 



1.8 X 10"^* m t = 9.6 X 10^ (10) 

Pap \ 100 days / ^ ' 

We therefore conclude that branching to^D dominates escape in the resonance lines up to at least 
~ 3000 days. 

A further important complication is that the resonance line photons at 1302 - 1306 A may 
be destroyed by photoelectric absorption by Si I, depending on the Si I abundance in the zone. If 
the probability of absorption dominates the branching probability, the 1302 - 1306 A photons are 
destroyed, and the branching io^D decreases drastically. This situation is similar to that earlier 
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discussed for helium. The probabiUty for this to happen can be estimated from equation (32) in 
Paper I, as 

P ~ A n ^Photo _ Ac. ^'^ 91 O'photo ngjl 

^photo ^ ^•'^ < 

Kline A'^ A21 92 noi 

= 1.2x10-3 f_^V" ^ ^1,) 
VIOOOK; noi ^ ' 

compared to the branching probabihty -Pbranch = 

3.0 X 10^6. Therefore, if nsii/riQi ^ 2.5 x 10"^^ 
photoabsorption dominates branching to the level. The triplet contribution to the level is 
therefore sensitive to the chemical composition in the oxygen-rich gas. Newer models of 15 and 
25 Mq stars by Woosley & Weaver (1995) give considerably lower silicon abundances in most of 
the oxygen zone. Further, the estimate in equation ([Til) is based on continuum absorption in the 
Doppler profile only. As we have discussed for helium, absorption also in the damping wings can, 
for complete distribution, increase the destruction probability above this value. This is confirmed 
by the model where we used the destruction probability from equation (33) in Paper I. In this 
model the non-thermal level decreased to a level of ~ 70% of that in Figure |lO[ As we have 
already remarked, however, we believe that the Doppler case is closer to the true situation of 
partial redistribution. 

We can now understand the origin of the somewhat paradoxical situation described above. 
The reason for the low contribution from the O - Si - S zone in the lOH model is the destruction 
of the 1302 A photons by the abundant Si I in this zone. If it had not been for this effect the 
contribution from the 3s state might have been similar in efficiency to that of the O - C zone, 
giving a flux proportional to the oxygen mass, and increasing the total flux by a factor ~ 3. Even 
in the O - Ne - Mg zone in the llEl model, the silicon abundance is relatively high, X(Si) = 
1.64 x 10"^, leading to the destruction of the 1302 A photons. In spite of the larger magnesium 
abundance, the absorption by Si I dominates the Mg I absorption, because the photoionization 
cross section for Si I is almost a factor of 100 higher than for Mg I. If for some reason we have 
over-estimated the continuum optical depth at 1302 A, e.g. because of too low ionization of Si I, 
the [O I] flux would increase. Charge transfer between O II and Si I could have this effect. In 
section we discuss this quantitatively, with negative result. 



To see the effect of a lower silicon abundance we show in Figure the [O I] light curve for a 
model where we have decreased the silicon abundance in the O - Ne - Mg zone in the llEl model 
to that in the O - C region, X(Si) = 1.3 x 10"^, compared to the original X{S\) = 1.64 x 10^^. 
We here see that the non-thermal part indeed increases by a factor of ~ 2. However, in spite of 
the smaller photoabsorption cross section of Mg I, in this model Mg I takes over the role of Si I, 
and most of the triplet contribution is also in this model quenched. An unwanted effect in this 
model is that there is a substantial increase in the flux between 600 - 850 days, which at this 
epoch destroys the agreement with the observations in the original model. 



The full extent of the photoabsorption can be seen in a model, also shown in Figure 12, where 
we have (artificially) decreased the photoabsorption of the A 1302 line to zero. Only in this model 
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do we get a non-thermal plateau close to the observed level. 

To estimate the dependence on the various parameters, and therefore more general models, it 
is of interest to consider a simplified model for the non-thermal [01] A A 6300, 6364 excitation. If 
we assume that the luminosity in these lines is determined only by the non-thermal excitation to 
the level, one can estimate the [O I] A A 6300, 6364 luminosity by the following argument. If 
J-y is the gamma-ray intensity, the absorbed energy in the oxygen component per unit volume is 
J-yK^po- Of this, a fraction e^xciD will give rise to direct excitations of the level. In addition, a 
fraction etrip will go into excitation of the level, via the triplet levels, as discussed above. The 
total energy fraction going into the level is therefore ei£) = e^^c^D + ^trip- Because collisional 
de-excitations at this epoch are unimportant, the total luminosity in the 6300, 6364 lines is given 
by 

Jeip K^{r) po{r)dV 

^6300,6364 - -f-j—lrVn TrUJV 

~ eiD — - — -^bol, (12) 

''"7, total 

where t^^q is the optical depth of the gamma-rays in the oxygen and r^_total is the total gamma-ray 
optical depth, both averaged over the ejecta according to the gamma-ray intensity. Both r^^o and 
'^7,totai depend on time as and equation (|l^) therefore shows that, except for the dependence 
of ei£) on the electron fraction, a fixed fraction of the bolometric luminosity is expected in the 
6300, 6364 lines, in qualitative agreement with the observations. 

From our Spencer-Fano calculations Cgxc.iD ~ 3.0 x 10"'^ (xe/10^^)~'^'^^^ for Xe in the range 
~ 3 X 10"^ - 3 X 10-2. At 1000 days we find Xe 7 x 10"^ so eexc,iD ~ 3.7 x 10"^ Our lOH model 
has Ty,o/''"7,totai ~ 0.18, so a fraction 1/6300,6364 /-^boi 6.7 x 10"^ of the bolometric luminosity 
should come out as A A 6300, 6364 emission. At 1000 days Lbd ?a 1.0 x lO^^ erg s'^ (Bouchet, 
Danziger & Lucy 1991), and -Z>6300,6364 ~ 3.8 x 10^^ erg s^^ (Danziger et al. 1991), so the observed 
ratio is -Z^6300,6364/^boi ~ 3.8 x 10"'^, while Menzies (1991) finds i6300,6364/^boi ~ (1.2 — 2.1) x 10^'^ 
at 1000 days. This shows that unless Xe <C 10^^ or Ty,o/'^7,totai ^ 0.18, direct excitation is not 
sufficient. 

The maximum contribution from the 35^5 level can be estimated by assuming that the 1302 
- 1306 A transitions have very large optical depths, and that photoabsorption is unimportant 
(see above). The triplet contribution to the luminosity in the 6300, 6364 lines is then 

T ~ \f ^ _912_ 1300 

-^6300, 6364, trip ~ [Jrccomb Cion ^onn ' ^exc,trip «orvnJ 

atot DoOO 6300 
^"•'^ Lboi. (13) 



"^7, total 



Here, ejon ~ 0.4 is the efficiency of ionization for O I (KF92), 03 the effective recombination 
rate to the "is^S level and atot the total recombination rate. From calculations by Julienne, 
Davies, & Oran (1974) we estimate that a^/atot ~ 0.22. The factor /recomb is the fraction of 
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the ionizations which recombine to O I. As we will see in section 5.2, /recomb can be very small 
if charge transfer of e.g. O II + Si I is efficient. Finally, eexc,trip is the fraction of the energy 
going into direct excitations of the triplet levels, of which most end up in the 3s^S state. In 
KF92 the total energy into excitations in the O - Si - S zone is ~ 0.1, of which ~ 42% is to the 
triplets, so eexc,trip ~ 4 X 10~^ at Xe ~ 10~^. The dependence on the electron fraction is roughly 
Cexc.trip oc x~^''^'^ , for 0.003 < Xg < 0.03. Therefore, the maximum efficiency for excitation from 
the triplets is etrip ~ 1.3 X 10-2 /recomb + 8.7 X lO^^. The maximum, total efficiency is then 

= fexc,iD 

+ etrip ~ 1.2 X 10-2(1 + 1.1 /.ecomb). With ''"7,o/''"7,totai ~ 0.18 and, depending on 
/recomb) wc find -i^6300,6364 ~ (2-5 — 4.9) X 10"^ -Z^bob in better agreement with the observations. 
This argument shows that the triplet contribution is sufficient, and most likely necessary, to 
explain the observed non-thermal level. 

The estimate above gives the luminosity for one particular hydrodynamical structure, with 
''"7,o/''"7,totai ~ 0.18. To check the sensitivity of this assumption we assume that all oxygen, 
with a total mass Mq, is located in a shell between velocities T^nin.o and Vmax,o- In this model 
T^{0) = Mo G/[4 vr (Kiax.o tf], where G = 3 [1 + (Kiin.o/Kiax.o) + (Klin,o/^max,o)^]"^ 
1 < G < 3. With this we get 

(14) 



-^^6300,6364 _ ei£, G Mq 



Lho\ 4 vr Vl^^ Q tl 

where t* is the time when the total optical depth, t^, total is unity. From light curve models 
Woosley, Pinto & Hartman (1989) find ~ 550 days. We therefore find that in the non-thermal 
phase 



^6300,6364 = 0.62 X 10-^ Q 



( ^'max,© \ ( ^ 



Mq 

1 Mq, 
-2 / J. \ -2 



V2000 kms-V V 550 days. 

(1 -I- 1.1 /recomb) -^^bol- (15) 

The main parameters determining the direct excitation of the lines are therefore the oxygen mass, 
the velocity interval of the oxygen mass, and to a weaker degree the electron fraction in the 
oxygen-rich gas. Of these the oxygen distribution is probably the main uncertainty, introducing 
an uncertainty of up to a factor of three. Equation (|l^) indicates that an oxygen mass of at least 
1 Mq is necessary to explain the observations (for G ~ 3). If most of the oxygen is in a narrow 
velocity range, G w 1, an oxygen mass of up to ~ 3 may be necessary. We caution, however, 
that due to uncertainties in both the atomic physics and the hydrodynamics this number is not 
firm. 

From this discussion it is obvious that there is a strong need in constraining the oxygen 
distribution from the observed line profile of the 6300, 6364 doublet. In Figure |l^ we show the 
A A 6300,6364 lines at 800 days for the llEl model (with our density distribution), including 
reddening and internal dust absorption. The observations are from Phillips et al. (1990). At this 
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epoch the lines are optically thin, and the doublet components can simply be added in a ratio of 
3:1. The ratio of the observed 6300/6364 lines is considerably smaller than the expected value of 
3.0 in the optically thin limit. A reason for this may be an additional contribution from the Fe I 
A 6361 line, which in our models is ~ 20% of the [O I] A 6364 line. The exact level of the Fe I 
A 6361 line is uncertain because of the uncertain UV-radiation field, and it is likely that it could 
explain the full discrepancy. 

In the velocity range 1000 — 2000 km there is a fairly good agreement with the observed 
slope of the line profile, indicating a realistic density distribution. At velocities ^ 2000 km s"^ 
the observations show a clear wing to ~ 3200 km s""^, considerably stronger than in the model. 
There may be two reasons for this. Either the oxygen abundance in the hydrogen envelope is 
larger than assumed in the model, or there is mixing of a small amount of processed oxygen to 
these velocities. The abundances in the hydrogen envelope is set by observations of the ring of SN 
1987A (Pransson et al. 1989, Sonneborn et al. 1996), and should be reliable. It would require very 
special conditions for the oxygen abundance in the envelope to be larger than that in the ring, 
which probably originates from a layer external to the envelope of the progenitor. We therefore 
believe that the most likely solution is the presence of some high velocity processed oxygen in the 
ejecta. 

The observed line profile is more peaked for V ^ 900 km s~^ than the model. Our innermost 
oxygen zone, containing ~ 0.6 M©, is at 770 km s^^, and it is clear that there is oxygen present 
down to at least ~ 400 km s~^. A better fit would be obtained if the same mass of oxygen-rich 
gas was mixed uniformly between 200 — 800 km s~^. The total flux should be the same and the 
total mass of the oxygen zone therefore similar to that in our model, ~ 1.9 M0, i.e. ~ 1.4 Mq of 
pure oxygen. 

In the same way as for the level, one expects non-thermal excitations to the level. In 
this case the contribution from the triplet levels is only ~ 2.5 X 10~3 of that to the level. The 

ratio of the non-thermal excitations to the and levels is £15 ~ 0.3 e^xc id &t Xe ^ 7 X 10~^ 
(KF92). The ratio of the A 5577 luminosity to the A A 6300,6364 luminosity is therefore 

L5577 .2972 1 

J = eis K-^^) ; 72972T- ^'-'^ 

i^6300,6364 5 5 77 ei^ + eis (f^) 

If excitations via the Ss^S level for the level are unimportant, ei£) = eexc,iD) ^^id the ratio 
of the [O I] A 5577 and [0 1] A A 6300,6364 lines is expected to be -^5577/^6800,6864 ~ 0.14. 
In the opposite case, when the triplet contribution is important, ei£) = egxc,iD + ^trip, and 
-^5577/-^6300,6364 ~ 2.5 X 10~^. The relative ratio can therefore give some information about the 
importance of excitation via the triplet lines of [0 1] A A 6300, 6364. The observed flux of the A 
5577 line is unfortunately uncertain because of blending with other lines. An approximate analysis 
of the spectrum at 804 days by Phillips et al. (1990) gives i^5577/i^6300,6364 ^ 0.3, as a strong upper 
limit. A more realistic limit, taking the blending by other lines (the 'continuum') into account, 
reduces this upper limit by a factor of two. 
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The O I A A 7774, 8446, and 9265 lines arise as a result of recombination, and in the case of 
the A 8446 line also by Bowen fluorescence by Ly/3 (Oliva 1993). The flux of the A 7774 line is 
on day 804 ~ 3 x 10~^^ erg cm~^ s~^. This is a factor of at least four stronger than observed in 
the day 804 spectrum by Phillips et al. (1990). Also the A 8446 line is over-produced by a similar 
factor. However, as suggested by Oliva, its absence can probably be explained by scattering by 
the Ca II IR triplet. The velocity difference to the closest component is only 1836 km s^^. The 
A 8446 photon may therefore contribute to the excitation of the Ca II triplet. We discuss this 
further in section 4.1.6. The A 9265 A line could be present in the spectrum. However, the peak 



of the observed feature is at 9234 A, which probably excludes it from being the O I line. The flux 
in this line is otherwise close to that expected for the [O I] line. The wavelength is consistent with 
Paschen 9-3 A 9234. The flux of this is, however, only expected to be 15% of that observed. 
There is also an excited multiplet in S I, As^S° — Ap^P, A A 9212 — 9237, coinciding with the peak 
of the line. This transition, which we do not include, can be expected to be strong because of 
recombination from the Si - S zone, and we think that this may be the most likely candidate for 
the line. Other strong quintet lines expected would then be the A 7696 and A 8694 lines. At 
the former wavelength there is a line with a flux of 40% of the A 9234 line, while the latter 
coincides with the Ca II triplet. 

Therefore, all the calculated O I recombination lines seem too strong, although the A 8446 
line could be scattered into the Ca II triplet. The K I resonance lines at 7664.9, 7699.0 A probably 
account for the P-Cygni line at ~ 7700 A. Its width is similar in velocity to the Na I D lines. It is, 
however, not likely that this doublet could scatter the A 7774 line, as this would probably result in 
a much stronger red peak of the line than observed. Instead, its total equivalent width is close to 
zero, as expected if it just scatters the 'normal' background emission. The wavelength difference 
is also too large. The effective recombination rates to the O I lines from Julienne, Davies, &: 
Oran (1974) are uncertain, but not by more than a factor of two, since the total recombination 
rate agrees within ~ 40% by that found by Chung, Lin, & Lee (1991). Instead, we think that 
charge transfer processes, not included in our model, are responsible for the quenching of the 
recombination lines. This is discussed in section ^.2| . 

Because of the large abundance and mass of oxygen, the bound-free emission continua of 

I may be observable, unless quenched by charge transfer. The strongest of these are the 
recombination continua to the 3p^P, 3d^D, and ScPDo levels, with edges at 4308 A, 8053 A, and 
8098 A, respectively. On day 800 the temperature in the oxygen zone is ~ 1500 K, corresponding 
to a width of ~ 660 A at 8000 A, and ~ 192 A at 4300 A. For the 3(fD and Sc^Dq emission we 
estimate a flux of 4.8 x 10~^^ erg s~^ cm~^, corrected for reddening and internal dust, at 8000 A. 
The observed continuum level at ~ 8000 A is ~ 9 x 10~^^ erg s~^ cm~^, consistent with the O 

1 continua. Also at 4300 A the model flux, ~ 1.5 x 10~^^ erg s~^ cm^^ is consistent with the 
observations, although blending with other lines makes it difficult to define the continuum level. 
Unfortunately, an unambiguous identification of these continua from the observations is difficult, 
although especially the continuum at 8000 A has no other obvious candidate strong enough. An 
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identification of this would make a direct determination of the O II fraction possible, and therefore 
to estimate the importance of charge transfer, as well as its temperature. 



4.1.5. Neon, Magnesium, Silicon, and Sulphur 

The neon abundance differs in the oxygen-rich zone by a factor of ~ 10 between the lOH and 
llEl models. This is a result of the different convective criteria employed in the two models. 
Consequently, the emission in the A 12.814 fim [Ne II] line differs considerably between the two 



models (Fig. 14). While the agreement with the lOH model is satisfying up to at least ~ 700 
days, the llEl model over-produces the line by a factor of ~ 2 — 3 at t 400 days. The last 
observations at ~ 730 days are considerably lower also in the lOH model. The uncertainty in the 
observation by Roche et al. (1993) at this epoch, however, is considerable. 



Mg I] A 4571 (Fig. ^5|) is completely dominated by the oxygen zone in both models. Because 
of the higher magnesium abundance in the llEl model, the flux is a factor of ~ 3 higher in this 
model compared to the lOH model. The agreement with observations is considerably better in the 
lOH model. The Mg I] line is interesting because it is dominated by recombination, rather than 
collisional excitation, as is usually the case. It is therefore not as sensitive to photoionization by 
the uncertain UV- field as one might think (see section |5.3| ) . The effective recombination rate of 
the line is, however, uncertain (see appendix in Paper I). Recombination dominates at all epochs 
over both thermal and non-thermal collisional excitation. 

Mg II A A 2795.5,2802.7 (Fig. |l|) is surprisingly weak, with a flux of ~ 10% of Ha, and 
similar to Mg I] A 4571. Before 750 days it is excited by thermal collisions, while at late time it is 
dominated by non-thermal excitation. The contribution is at early time dominated by hydrogen 
and helium regions, while in the non-thermal phase the oxygen-rich region dominates. Because of 
a lower Mg II fraction at late time, the line is in this phase weaker in the llEl model compared to 
the lOH model, despite a higher magnesium abundance in the former. Unfortunately, resonance 
scattering by the many metal lines in the UV makes the observed luminosity at ^ 1200 days of 
this line uncertain. The same is true for the Mg I A 2852 resonance line. In section ^ we compare 
it to the HST observations at 1862 days. 

The 1.64 /xm line is most likely a blend of the [Si I] A 1.6454 ^um and [Fe II] A 1.6435 //m 
lines. The observed, relative contributions are, however, uncertain. The [Si I] A 1.6073 //m line 
should be a factor 2.84 weaker than the [Si I] A 1.6454 //m line. Also this line is blended with 



[Fe II] and Br 13. Figure 17 shows the calculated fluxes in the 1.64 //m feature, together with the 
individual contributions. During most of the evolution the total flux is over-estimated by a factor 
2-3. The reason for the early over-production may be the large contribution from the O - Si ~ S 
region, which dominates at that time. A lower abundance of silicon in this region, may bring the 
early part into agreement with the observations. At t ^ 1000 days the [Fe II] contribution from 
the hydrogen region dominates in all models. The [Si I] A 1.0991 /xm line is blended with Pa7, 
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and also with the wing of the He I and [S I] hnes, discussed below. Therefore we do not discuss it 
in any detail, and just note that the fluxes given by Meikle et al. (1993) for this hne are somewhat 
lower than those from the lOH model. The Pa7 flux should be small compared to the [Si I] line 
(Fig. Q), but the contributions from the He I and [S I] lines large. 



As we have already discussed in section 4.1.2, the strong 1.08 fim feature is most likely a blend 



of Pa7, He I A 1.0830 fim and [S I] A 1.0820 fim. At t ^ 400 days there is too high a luminosity 
in this feature in the lOH model (Fig. ^). The good agreement with the He I A 2.058 fim light 
curve gives us confidence that the He 1 A 1.0830 /im is fairly accurate. The [S I] line is therefore 
likely to be too strong in the lOH model. This is mainly a result of the high sulphur abundance in 
the O - Si - S region in the lOH model. 

[S I] - 3p^ at 1.0820 ^m and [S I] 3p^ ^Pi - 3p^ at 1.1306 /im originate from 

the same upper level. The probability for emitting a A 1.0820 /im photon is 77 %, and 23 % for a 
A 1.1306 /im photon. From observations by Meikle et al. (1993) we estimate the flux in the feature 
at 1.13 /im to ~ 10~^^ erg s~^ cm~^ at day 695. Another possible candidate to this feature is O I 



A 1.1287 /im. As was mentioned in section 4.1.4, this line may arise as a result of the 1- Ly/3 



fluorescence mechanism (Oliva 1993). We have not included this process in our calculations, and 
can therefore not estimate the efficiency of this mechanism. Assuming the 1.13 /im feature to be 
due to [S I] A 1.1306/im results in a maximum estimated flux of the [S I] A 1.0820 /tm line of 
~ 3 X 10~^^ erg s^^ cm~^. In our calculations the flux in [S I] A 1.0820 fim is 8 X 10-^2 and 
~ 3 X 10~^^ erg s~^ cm~^ for the lOH and llEl models, respectively. Also from this line the llEl 



model is therefore favored by the [S I] emission. As mentioned in section 4.1.4 , there may be other 
allowed lines from S I in the spectrum, not included in our model. 



4.1.6. Calcium 

The light curves of the [Ca II] A A 7291, 7324 lines and the IR-triplet are shown in Figures 



18| and |19|, together with contributions from the different components. The overall agreement of 
the A A 7291, 7324 light curve with observations is for the lOH model quite satisfying, while the 
lines are severely under-produced in the llEl model. At t 900 days there is a clear tendency to 
under-produce the luminosity in both models. This is still more pronounced with the triplet, but 
here the contributions from the [C I] A 8729 and O I A 8446 lines also have to be included. 

An interesting point is the origin of the emission at various phases, which is sensitive to the 
explosion model and differs between the lOH and llEl models. At no epoch in either model does 
the Si - S component, where the explosively synthesized calcium resides, dominate. Instead, either 
the oxygen region or the hydrogen-rich gas contribute most of the emission in the lines. 

We agree with Li &: McCray (1993) that little of the Ca II emission comes from the Si - 
S zone, although our conclusion about the origin differs somewhat from that of Li & McCray. 
These authors find that the Ca II lines all originate from primordial calcium in the hydrogen and 
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helium-rich matter, with a filhng factor of ~ 0.1. Based on the hkely calcium mass, they argue that 
the newly synthesized calcium captures too small a fraction of the total gamma-ray flux to explain 
the observed line strengths. Their argument rests on the assumption that calcium is the dominant 
element in the newly synthesized gas. This is not the case, and nucleosynthesis models show that 
the calcium fraction is only ~ 0.05, with most of the mass in silicon and sulphur. Because of the 
efficiency of Ca II as a coolant (e.g., Fransson &: Chevalier 1989), most of the cooling of this mass 
may in fact be done by Ca II lines. The reason for the dominance of the other regions is therefore 
the low total mass of the Si - S zone, compared to the other regions. In the lOH model the mass 
of this region is ~ 0.30 Mq, and a fraction of only 5 x 10~^ of the gamma-rays are absorbed in 
this gas. 

As found by Fransson & Chevalier (1989), Ca II may be a strong coolant of the oxygen-rich 
gas, if mixing of calcium during the hydrostatic pre-explosion burning is efficient. Crudely, the 
critical calcium fraction for Ca II to dominate the cooling is X{Ca) ~ 10~^. The extent of this 
mixing is determined by the convection criterion used, as well as the extent of over-shooting. In 
the lOH model, where the Ledoux criterion was used and over-shooting included, this mixing 
is important, resulting in an abundance of ~ 6 x 10~^ in the O - Si - S zone, by number. 
Consequently, at t ^ 400 days the O - Si - S region dominates the A A 7291, 7324 contribution, 
while at later time the hydrogen component within the core, with unprocessed calcium, dominates. 
Already at 700 days the temperature in the hydrogen-rich gas is more than a factor of two higher 
than in the O - Si - S region, explaining the higher contribution to the Ca II emission from the 
former component. The slow adiabatic decrease of the hydrogen temperature makes the decline 
of the Ca II lines from this component relatively slow. The temperature effect is seen even more 
clearly in the triplet emission, which is dominated by the hydrogen core component at all times. 
The llEl model, on the other hand, has very little calcium mixed with the oxygen. 

Li &: McCray find that in order for the Ca II lines to be at the observed level after ~ 400 
days, radiative pumping of the H and K lines is necessary. As is seen in Figures ^ and 19, our 
light curves agree fairly well with observations up to ~ 800 days, if all lines coinciding with the 
7300 A lines and IR-triplet are included. The O I A 8446 line has a velocity difference of 1836 
km from the Ca II A 8498 component, and 3372 km s^^ from the stronger A 8542 line. 
These lines are optically thick in both the hydrogen and oxygen components, as well as in the Si - 
S - Ca region. In the H-core the A 8498 component is optically thick up ~ 750 days, while the 
A 8542 line is optically thick in the H-envelope at ~ 3400 km s~^ up to ~ 350 days. The O I 
A 8446 photons will therefore be scattered and emerge as Ca II triplet emission up to ~ 800 days. 
This justifies the inclusion of the line in the flux from the Ca II triplet, although in section 5.2 we 
argue that charge transfer probably decreases the 8446 flux dramatically. 

At times later than ~ 800 days there is a clear deficiency. This is especially clear when 
observations from CTIO at 953 - 1149 days are included (Suntzeff et al. 1991). The observed flux 
in the 7300 line on day 1046 is a factor of ~ 10 higher than in our models, and the 8600 complex a 
factor of two, including O I A 8446, which may be questionable. Our models do include radiative 
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excitation, but only by the continuum. On day 1046 the two-photon continuum flux within the 
region covered by the H and K lines, ~ 170 A, is only ~ 6 x 10~^^ erg s^^ cm^^, which is much lower 
than the total, dereddened flux in the 7300 A lines and IR-triplet, ~ 2.9 x 10~^^ erg cm~^. 
If we, however, for the same day extrapolate the observed average, dereddened flux on either side 
of the H and K lines, ~ 1.3 x 10~^^ erg s~^ cm~^ A , to the wavelength of these lines, we can 
estimate the total absorbed flux by the H and K lines to ^ 2.3 x 10^^^ erg s~^ cm~^, close to the 
total flux in the 7300 A lines and IR-triplet. Although the flux in the interior of the ejecta may 
differ from the escaping flux, this strongly suggests that radiative excitation in the H and K lines 
is responsible for the flux in the 7300 A lines and IR-triplet at late epochs. The fact that the 
fluxes of these lines follow the bolometric flux at a level (-L7300 + -^iR)/-^^boi ~ 7 x 10"'^ does not 
mean that non-thermal excitation is necessary, as for the [O I] lines. If the U-band flux follows the 
bolometric, as is a fair approximation to the observations after ~ 800 days (Suntzeff et al. 1991), 
radiative excitation will also result in a nearly constant (-L7300 + -^^iR)/-^boi- In fact, using the 
efficiency for non-thermal excitation in KF92, one can estimate the non-thermal contribution to 
^ 10~^ (Mgi_g/0.3 Mq) of the bolometric, much less than observed. 

The relative intensities of the Ca II H and K lines, the A A 7291, 7324 lines and the IR-triplet 
are discussed under various conditions by Ferland & Persson (1989), Pransson &: Chevalier (1989), 
and Li & McCray (1993). These authors find that above Tie ~ 10^ cm~^ the H and K, as well as 
the IR triplet, are in LTE, and therefore mainly depend on temperature. At temperatures below 
~ 4500 K both the H and K and the IR triplet decrease rapidly with temperature. At times when 
radiative excitation dominates, and if thermalization of the lines can be ignored, the ratio of the 
7300 A lines and IR-triplet is expected to be ~ 8600/7300 = 1.18, close to the observed value of 
~ 1.4 at 1046 days. 



4.1.7. Iron, Cobalt, and Nickel 

In our model Fe I- IV are treated as multilevel atoms with 121 levels for Fel, 191 levels for 
Fell, 110 levels for Felll, and 43 levels for FelV. Fe V is included only with its ground state. 
The atomic data, which includes new IRON project data, are discussed in the appendix in Paper 
I. The total recombination coefficients (Shull & Van Steenberg, 1982), and also the fractions of 
the radiative recombination rates going to the ground states (Woods, Shull, & Sarazin, 1981), are 
probably reasonably well known. However, a major uncertainty is the individual recombination 
rates to the excited levels of Fel, and Fell, which are largely unknown. We will discuss the 
consequences of this later. 

In Figure |2^ we show the light curves of the most important Fe II lines. The observational 
data are taken from Erickson et al. (1988), McGregor (1988), Meikle et al. (1989), Moseley et 
al. (1989), Haas et al. (1990), Varani et al. (1990), Spyromilio et al. (1991), Dwek et al. (1992), 
Jennings et al. (1993), Colgan et al. (1994), and Bautista et al. (1995). 
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In general we find agreement with observations to within a factor of two, or better. It 
should be noted that the agreement with observations at epochs earlier than ~ 400 days is very 
good, with a possible exception of the heavily blended 1.53 /xm line. At later epochs most of the 
calculated lines without correction for dust absorption are a factor ~ 4 stronger than observed. 
Including this correction with a factor ~ 40%, as found by Lucy et al. (1991), gives a greatly 
improved agreement, although they are still up to a factor of ~ 2 over-luminous. The importance 
of dust was earlier noted by Colgan et al. (1994) for the A 25.99 /im line at 640 days. 

The contributions to especially the A 7155 line and the near-IR A 1.26 /^m and A 1.53 /xm 
lines are interesting. From being dominated by the iron core at t ^ 600 days, the contribution 
from unprocessed iron in the hydrogen component becomes the most important after this epoch. 
The same is true for the far-IR A 17.94 /im and A 25.99 ^m. lines, although the transition for these 
occur later. For these lines the hydrogen component is important throughout the whole period, 
and later than ~ 800 days dominates the iron core contribution. The rapid drop in the different 
lines from the core can clearly be seen in Figure 6 of Paper I. The dominant hydrogen contribution 
to the A 25.99 /xm line at late time is in line with the KAO-observation at 1153 days by Dwek et 
al. (1992), who find an emitting iron mass of only 3.8 x 10~^ Mq, consistent with primordial. 

An estimate of the luminosity of the A 25.99 /xm line coming from unprocessed gas in the 
hydrogen to oxygen regions is given by 



— 2.7 X 10'^^ ^rr,L 



X(Fe) \ f M 



1.7 X 10- V Vl Mq 

g-554 K/T 



where M is the total mass of hydrogen, helium or oxygen-rich gas. The optical depth of the 
A 25.99 ^m. line in these components is 



2000 kms-V V4000 Ky V500days 



(18) 



where V and T are the average velocities and temperatures, respectively, of the unprocessed gas, 
and fi is the filling factor of the component. We assume that T ^ 554 K. At 400 - 600 days 
the typical temperature in the hydrogen and helium-rich gas is ~ 4000 K (Paper I). Further, 
our line profile fits show that most of the hydrogen and helium-rich gas have velocities within 
~ 4000 km s~^. Therefore, unless clumping is high, the A 25.99 fim line from the unprocessed 
iron is likely to be optically thin, as is also argued on observational grounds by Haas et al. (1990). 
The total luminosity of the A 25.99 /im line on day 407 was ~ 1.3 x 10^'' erg s~^, decreasing to 
~ 6.3 X 10^6 erg s"^ on day 640 (Haas et al. 1990, Colgan et al. 1994). With M ;^ 10 Mq, it is 
clear that a large fraction of the A 25.99 /im line may originate from unprocessed iron. As Figure 



20 shows, the same applies to the other [Fe II] lines, including the A 1.26 /im and A 17.94 /im line. 
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The disappearance of the iron core contribution to the near-IR [Fe II] hnes was first noted by 
Spyromiho & Graham (1992), who correctly attributed this to the IR-catastrophe. From their last 
observation on day 734 they infer an Fe II mass of ^ 0.01 Mq, and attribute the rest to iron too 
cold to emit in the near-IR. 

While our light curves agree well with those of Li, McCray Sz Sunyaev (1993) for t ^ 700 
days, a major difference between our models and those of Li et al. is in the behavior of the Fe 
II light curves at late time. The models by Li et al. fail to explain the emission from the Fe II, 
Co II, and Ni II lines for t ^ 2 yr, with too rapid a drop of especially the optical and near-IR 
lines at ~ 600 days. As a solution they suggest that photoionization by the UV continuum, 
especially two-photon emission, from helium surrounding the iron may solve the problem. We 
do not experience this problem, because of the contribution to these lines from unprocessed iron 
in the helium and hydrogen region. In addition, we do include photoionization in the core, and 
as we discuss below, this may actually give too large an effect on the state of ionization. If 
photoionization really was needed, there may be a more local source of UV-photons; according to 
the lOH and llEl models there is a substantial amount of helium mixed microscopically with the 
iron, roughly equal abundance by number (see Tables 1 and 2 in Paper I). This could result in a 
substantial local UV-radiation field in the iron clumps (see below). 

The hght curves in Figure |l] of [Co II] A 10.52 /im, and [Ni II] A A 6.634, 10.68 ^um all 
agree as well as can be expected, given the quality of the atomic data. Up to ~ 700 days the 
contribution from the Fe - He core dominates all lines. After this epoch the Si - S contributes 
most of the flux of [Ni II] A 6.634 //m up to ~ 1100 days. Primordial nickel dominates [Ni II] 
A 10.68 /xm already after ~ 600 days, similar to the [Fe II] lines. The [Co II] line drops rapidly 
after ~ 900 days because of the low primordial cobalt abundance. 

While agreement with the Fe II, Co II, and Ni II lines is satisfying, the Fe I lines are weaker 
in these models by several orders of magnitude. This is seen in especially the [Fe I] A 1.44 /xm 
emission. The under-production of the [Fe I] A 1.44 ^um line is an ionization effect, resulting from 
a very low predicted Fe I abundance in all regions. The typical Fe I fraction in the iron core at 
400 - 600 days is only (1 — 4) x 10~^, with most of the ionization resulting from photoionization 
by Fe II recombination lines and the He I two-photon continuum. A similar effect is seen in the 



[Ni I] A 3.119 /im line, which is under-produced by a similar magnitude. In section 5^ we show 
that changes in the UV-field can have dramatic effects on the fluxes of these lines. 



5. UNCERTAINTIES IN THE CALCULATIONS 
5.1. Filling Factors 



The filling factors we employ for the various components are uncertain. To test the sensitivity 
of our results to the assumed filling factors, we have run a set of models where we have varied all 
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filling factors by a factor of two in either direction (always assuring that the total is one). The 
range we have investigated is therefore /h = 0.075 — 0.30, fue = 0.13 — 0.60, /o-c = 0.03 — 0.12, 
/o-Si-^s = 0.06 - 0.24, /si„s = 0.0085 - 0.034, and /pc = 0.17 - 0.70. 

With regard to the temperature evolution, we find a relatively small change up to the time of 
the IR-instability. The epoch when this sets in varies in the metal-rich regions by ~ 100 days, over 
the whole range of filling factors. The IR-instability occurs earlier when the filling factor increases, 
i.e., the density decreases. This affects the epoch of the steepest part in the light curve of the 
different lines. Because of the importance of adiabatic cooling, the temperature of the hydrogen 
and helium components in the core are hardly affected at all. 

Most of the hydrogen lines are insensitive to the filling factor. An exception is H/3, which 
in the period 200 ~ 400 days increases by a factor of 2 when /h decreases from 0.3 to 0.075. 
The insensitivity of the hydrogen lines to the filling factor is in some disagreement with Xu et al. 
(1992). 

The [0 1] A A 6300, 6364 lines show little dependence on the filling factors in either the 
oxygen or hydrogen and helium-rich regions. The non-thermal part is also unaffected by this. 
The discrepancies in the light curve of this line are therefore not likely to be connected to the 
assumed density, but rather to the chemical composition or charge transfer effects. Also [Ca II] 
A A 7291,7324 show weak dependence on the filling factors. The exception is the contribution 
from the hydrogen core, which in the period 400 - 800 days decreases by a factor of up to three 
between /h = 0.075 and /h = 0.3. The lower filling factor would give a better representation of 
the light curve. The contribution from the O ~ Si - S zone, however, differs considerably less as 
/o varies. Roughly the same is true for the Ca II triplet lines. 

While the optical and near- IR Fe II lines, e.g., A A 7155 A, 1.26, 1.53, 1.64 //m, are insensitive 
to both fpe and /h, the far-IR lines show a higher sensitivity. Of the individual contributions to 
the A 17.94 ^m and A 25.99 /im lines the hydrogen core contribution is only weakly dependent 
on /h. The Fe - He contribution, however, decreases by a factor nearly proportional to /pe- These 
conclusions agree with equation ( p^ and equation (43) in Paper I. Contrary to Li, McCray, & 
Sunyaev (1993) we find that the observations agree better with a fairly low filling factor, /pe ~ 0.2. 
The reason may be our inclusion of regions other than the iron-core, and the additional Fe II 
emission from these. 

Summarizing this discussion, we find that plausible variations of the filling factors have 
rather small effects, at the factor of two level, or less. Except for perhaps the far-IR [Fe II] lines, 
it is likely that other uncertainties, e.g., the hydrodynamics or abundances within the different 
components, are at the same level or worse. We are therefore somewhat cautious of drawing any 
far-reaching conclusions about the filling factors from this. 
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5.2. Charge Transfer 

Charge transfer is important for both the ionization balance and the hne emission. In 
particular, lines and continua arising as a result of recombination can be severely affected. Earlier 
we noted a marked over-production of the O I recombination lines. As was shown in Paper I, the 
O II fraction is, however, sensitive to charge transfer with Si I, with an uncertain rate. To see the 
effect of this in more detail we have varied the rate of Si I + O II — > Si II + O I in the region 
10"^'^ — 10~^ cm^ s~^ in our models. We find that the O I A A 7774, 8446 and 9265 recombination 
lines, as well as the O I continua, all decrease by factors of 10 - 100, corresponding to fluxes well 
below the upper limits from the spectra. This is obvious from the O II curve in Figure 10 of Paper 
I. Charge transfer between O II and Si I would therefore solve the problem with the too large O I 
recombination fluxes in the previous models. At the same time the problem with the non-thermal, 
flat part of the [0 1] A A 6300, 6364 light curve is exaggerated, because the contribution from the 
triplet levels, which to a substantial part are fed by recombination, decreases too. 

Because of the small change in the Si I fraction, [Si I] A 1.6454 //m does not change 
appreciably. In Paper I we noted that the increase in the Si I + O II charge transfer lead to a 
decrease in the Mg I fraction. This does, however, not change Mg I] A 4571 substantially, because 
of the dominance of recombination. Because of the large contribution to the 8600 A feature from 
the O I A 8446 line without O II + Si I charge transfer (Fig. |l^) , the flux of this decreases in this 
model substantially, later than 800 days. Photoexcitation can probably compensate for this. 

Another interesting consequence of the O II + Si I charge transfer is that the Fe I fraction 
in the O - Si - S zone increases by several orders of magnitude, although still too low to give 
Fe I lines of sufficient strength. The reason for this is that the UV flux in this zone to a large 
extent is determined by O I and Si I recombinations. Both decrease dramatically, while the Mg I 
recombination radiation is not energetic enough to ionize Fe I. In principle, other charge transfer 
processes in the other zones could possibly increase also the total Fe I flux. 

We have also tested the influence of charge transfer from excited states of He I, discussed by 
Swartz (1994). However, here we find that they only contribute a few percent of the total He I 
ionization. Charge transfer from the ground state is also unimportant for helium. 

5.3. Photoionization 

Resonance scattering in the hydrogen dominated envelope has been discussed by Li & 
McCray (1996), using a Monte-Carlo model. Although their results are sensitive to their assumed 
temperature, the results may be indicative to the order of magnitude in the intensity. As input 
spectrum they use a pure He I two-photon continuum, but do not include other sources like the 
H I two-photon continuum or line emission, like Lya, 1 A A 1356, 1641, and Mg II, as well 
as recombination emission from Si I and Mg I. Nevertheless, these calculations illustrate the 
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importance of the UV scattering. 

In Li & McCray's model at 200 days the emergent flux is down by a factor of ~ 10, while the 
lower temperature and density at 800 days give a factor of ^ 3 suppression. The intensity inside 
the ejecta should be larger than the emergent intensity. On the other hand, scattering by the 
iron-rich material in the Fc - He bubble, not considered by Li & McCray, may further increase 
the scattering, and decrease the intensity. Resonance scattering will decrease with time, as is 
shown by the increasing UV flux in the lUE-band of SN 1987A (Pun et al. 1995), as weU as the 
emergence of clear lines in the HST-spectra at ^ 1800 days (Wang et al. 1996, Chugai et al. 1997). 

In addition to resonance scattering, which mainly redistributes the emission in wavelength, 
dust absorption may have strong effects for the UV intensity. The resonance trapping in the UV 
increases the path length, compared to the optical range, and with a covering factor of ~ 40% the 
probability for absorption may be very large. Experiments with a simplified Monte Carlo model, 
similar to that in Pransson (1994), confirms this. 

In KF92 the different contributions to the UV cmissivity were discussed in detail, using 
arguments based on energy input, optical depths and atomic physics. Most of the results in KF92 
are confirmed by our calculations, although quantitatively the emissivity difi'ers from KF92. The 
most important reason for this is the different mass distribution in the models here, based on 
observed line profiles, and in KF92. 

In the hydrogen zone most emission comes from H I, with Lya and two-photon emission being 
most important. When the optical depth in Lya is sufficiently large the photons scatter many 
times, and finally emerge in the two-photon transition from the 2 state to the ground state 
(Xu et al. 1992, KF92, Fransson 1994). In the hydrogen-rich regions in the core we find that the 
two-photon continuum dominates over Lya (i.e. A2-f > ALya/3csc) at all times between 200 - 2000 
days. In the hydrogen envelope, where the density is decreasing outward, Lya dominates in the 
outer zones at 500 days, and by 1200 days the Lya emission is dominating in the entire envelope. 
At 500 days 45 % of the total energy absorbed by the hydrogen envelope is emitted in the H I 
two-photon continuum. 

In KF92 it was found that ~ 50% of the energy absorbed by helium is re-emitted as 
He I two-photon emission. Of the emission coming out above 11.26 eV, most is absorbed by 
photoelectric absorption by the C I in the helium zone. This is caused by the large carbon 
abundance in this zone. A large fraction of this is in turn re-emitted as C I recombination 
lines in the UV. Altogether, it was estimated that ^ 40% of the energy was emitted in the UV 
between 1100 - 3646 A. This high efficiency coupled to the large relative fraction of the energy 
deposited in the helium zone, means that the helium zone can be the most important source for 
the UV-radiation field. 

In these calculations we find that compared to KF92 a considerably smaller fraction of the 
deposited energy goes into the helium zones, only ~ 4 — 8%, (compared to ~ 30% in KF92). This 
is mainly due to our redistribution of the helium to higher velocities. Of the total emission from 
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the helium zones, ~ 30 - 60% is emitted in the UV (A < 3646 A). Of this, ~ 40% comes from the 
He I two-photon continuum after ~ 500 days. The 1100 - 3646 A UV-efficiency is similar to the 
number above, ~ 37% of the absorbed energy. As in KF92, we find that the 584 A line is mainly 
recycled into the 2.058 /um line, plus two-photon continuum, because of the higher probability of 
branching into the 2.058 fim line, compared to escape or continuum destruction in the 584 A line. 
The He I and H I two-photon continuum photons constitute the dominant sources for the UV-field 
in the ejecta. Up to ~ 600 days other strong UV-lines are Mg II A 2800, C II] A 2326, and C I 
A A 2966, 2968. 

In the oxygen region recombination emission to the ground state of Si I (~ 14% of the 
absorbed energy) and Mg I (~ 3%) at 1521 and 1621 A, respectively, give most of the emission. 
Also recombination to the excited states in Si I is important (~ 9%). Because our recombination 
model of Si I is crude, we can not predict the detailed emissivities of these lines, although most 
should emerge at ~ 2000 A. 

The UV fraction from the Fe - He rGgion is coniparcitivGly low, L(^X <C 3646 A.)/-Lf;otai 

« 0.14. 

Most of this emerges as He I recombination radiation to excited states, and He I two-photon 
continuum emission (~ 2.6%). Because of the large optical depths most Fe II emission is degraded 
into optical lines. 

To investigate the uncertainty in the UV radiation field we have made one calculation where 
all photoionization processes from the ground states of the low ionization potential elements, Na 
I, Mg I, Si I, Ca I, Fe I, Ni I and Co I, are switched off. 

Without photoionization we now get an Fe I fraction of ~ 17% at 500 days, while with 
photoionization the Fe I fraction is ~ 10~^. This illustrates the sensitivity of the neutral fraction 
of low ionization potential elements. The Co I and Ni I fractions increase correspondingly. 

Our [Fe I] and [Ni I] light curves (Fig. [2^ ) are without photoionization only a factor of two 
lower than the observations, but given the quality of the atomic data and the sensitivity to both 
UV-field and filling factor, it is not too surprising. Because the Fe I fraction is sensitive to the 
assumed filling factor of iron, with n(Fe I)/n(Fe) oc fp^ for n(Fe I)/n(Fe) <S 1, a smaller filling 
factor could increase the fluxes. The [Fe II] lines decrease by ~ 25%, corresponding to the decrease 
in the Fe II fraction, improving the agreement with the observations somewhat. 

The Mg I] A 4571 line decreases by a factor of 2 - 3 for t ^ 700 days, in better agreement 
with the observations. The reason is that radiative recombination is responsible for ~ 98% 
of the flux of this line at 800 days. Because the fraction of Mg II at 800 days has decreased 
from X{Mg II) « 1.00 to X{Mg II) ^ 0.28, the Mg I] flux decreases. At t ^ 500 days the 
Mg I] flux is instead higher in the model without photoionization. This is consistent with a 
large recombination contribution to the line, which does not vary much in the two models, 
because X(Mg II) ~ 1.0. In the non-photoionization model X(Mg I) ~ 1.1 x 10^^, while it is 
~ 7.3 X 10~^ with photoionization. Without photoionization there is therefore a substantial 
collisional contribution to the line, explaining the increase. 
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Summarizing this, there are strong indications that the UV-field, below at least ~ 1600 A, in 
the ejecta is substantially weaker than in our standard models, because of resonance scattering, 
in combination with branching into the optical, or, especially efficient, dust absorption. We 
emphasize that this uncertainty mainly affects ions with low ionization potentials, and abundances. 
Dominating ions, like H I, He I, O I, Mg II, Ca II and Fe II, are only affected marginally. 

6. DISCUSSION 

Summarizing these calculations, we find both successes and failures. Starting with the former, 
we find a general good qualitative agreement with the observed evolution of most lines. 

The general behavior of the hydrogen and helium lines is well reproduced, although this has 
already been shown by previous calculations (Xu et al. 1992, KF92). Here, we have added the 
time dependence, which is necessary for t ^ 800 days. In KF92 we found a deficit in the flux at 
late time. This has now disappeared, partly because of the freeze-out (see also Chugai et al. 1997). 
In addition, we include a more realistic temperature dependence, as well as a more consistent 
calculation of the intensity of the radiation in the Balmer continuum, which in turn determines 
the Balmer line fluxes. 

The effects of including the time dependence are most pronounced for lines excited by thermal 
collisions and recombination lines, coming from regions where adiabatic cooling and freeze-out of 
the ionization are important. These effects cause the electron density, as well the temperature, to 
be larger than in steady state (Paper I). Lines dominated by recombination are affected mainly by 
the higher electron density compared to steady state, increasing their luminosities. Lines formed 
mainly by non-thermal excitation arc almost unaffected. Consequently, the time dependence is 
seen clearest for the [Fe II] lines at late stage, when the hydrogen and helium regions dominate. 
In particular, the optical and near IR-lines, e.g. A A 7155 A and 1.26 //m, increase at 1000 
days by factors 3-4, compared to steady state. Also the A A 17.94, 25.99 /xm lines increase by 
50 — 100%. The H I and He I lines, however, increase by only ^ 50%, at the same epoch, while [O 
I] A A 6300, 6364 is almost unaffected. 

A new feature of this paper is the use of the line profiles for constraining the distribution 

of the various components. This was discussed already in Fransson & Chevalier (1989) and in 
FK92, but mainly from a qualitative point of view. Here, wc have shown how both the hydrogen, 
helium and oxygen distributions can be determined from the observations. In particular, the 
penetration of a substantial hydrogen mass to ^ 700 km has been established. This has earlier 
been proposed in connection to light curve calculations (e.g., Shigeyama & Nomoto 1990). Also 
oxygen and helium are mixed to low velocities, in the case of oxygen to ^ 400 km s~^. In the 
other direction, synthesized iron is mixed outwards to at least 2000 km s~^. Conclusions about 
mixing of newly synthesized iron to even higher velocities should, however, be taken with caution, 
because of the strong contribution of primordial iron to the lines. Because of the low resolution of 
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the spectra used, it would be desirable to repeat this analysis with data of higher quality. 

When discussing the masses of individual elements, the luminosity of the different lines are 
of obvious interest. However, the deposition of the gamma-rays depends on the distance from the 
radioactive source. While the deposition per unit mass is fairly uniform inside the radioactive 
source, i.e. inside ~ 2000 km s~^, it decreases like outside. The mass outside the radioactive 
source is therefore not well constrained, unless the wings of the lines are well reproduced. Neglect 
of this leads to an under-estimate of the mass. 

Of the masses determined, the hydrogen zone mass is probably the most accurate, both 
because of good observations, and because of well-understood physics. The total mass we find, 
~ 7.7 Mq, in the hydrogen component is probably accurate to ±2 Mq. Of this ~ 3.9 is pure 
hydrogen. 

The quality of the line profile of the least blended of the He I lines, A 2.058 /xm, is 
unfortunately not ideal. Especially the wings are not very well defined. Therefore, although the 
mass of the helium-rich component inside ~ 4000 km s~^ is consistent with ~ 2 Mq, there could 
be a substantial additional mass at velocities ^ 3000 km s^^. There is also an uncertainty in the 
helium mass connected to the form of the continuum destruction probability used. As illustrated 
in Figure ^, there may, depending on this, be a factor of up to five uncertainty in the flux at early 
time. At ^ 500 days this uncertainty is considerably smaller, and this part of the light curve is 
therefore better suited for helium mass determinations. We again note that ~ 3.9 Mq of helium 
resides in the hydrogen-rich gas. 

In both the llEl and lOH models most of the [C I] emission originates from regions outside 
of the O ~ C shell. Although this is not the case in the model with the He - C shell replaced by 
a He - N shell, the influence of CO as a coolant in the O - C shell makes mass estimates based 
on the [C I] lines questionable. In the model where we replaced carbon by nitrogen in the helium 
shell, as we argue is indicated, the [N I] A 1.04 /im line is dominated by the He - N region. The 
nitrogen mass in this model is ~ 3.4 x 10^^ Mq, of which ~ 2.4 x 10~^ Mq is in the He - N zone. 

The strength of [Ne II] A 12.814 //m is in the lOH model dominated by the helium zone, 
while in the llEl model the O - Ne - Mg zone takes over. The flux in the [Ne II] line was in the 
latter over-produced by a factor of 2 ~ 3 later than 500 days . Therefore, the neon mass in the 
llEl model, ^ 0.25 Mq, is too large, and can be seen as a solid upper limit. The neon mass in the 
oxygen zone of the lOH model, ~ 2 x 10^^ Mq, together with ~ 4 x 10~^ Mq in the helium zone, 
gives an acceptable fit to the light curve. Allowing for a factor of two uncertainty in the model, 
the neon mass in either of these regions, but not both, can possibly be increased by a factor of 
two. A best estimate is therefore M(Ne) f» 6 x 10^^ Mq, and most likely ^0.1 Mq. 

Masses based on Mg I] A 4571 have earlier been questioned, based on the sensitivity on 
the far-UV radiation (c.f., Fransson & Chevalier 1989). The fact that the line is dominated by 
recombinations from Mg II makes this determination considerably more reliable. The magnesium 
mass is likely to be similar to that in the lOH model, ~ 2.2 x 10~^ Mq. The main uncertainty in 
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this comes from the adopted effective recombination rate of the A 4571 hne. 

As we have discussed, and agreeing with Li & McCray (1993), calcium is not well constrained, 
because of the dominance of the zones outside of the Si - S zone for the emission in the Ca II lines. 

The iron mass should be close to the value used here, certainly within a factor of two. This 
mass, as well as the ^^Co mass, is however, more accurately determined from the bolometric light 
curve. The amount of stable nickel, ^^Ni and ^'^Ni, is in our model ~ 0.006 and given the 
uncertainty in the atomic data for Ni II, this is consistent with the observations. 

We now turn to a number of more specific points. Of special interest is to examine the 
sensitivity of the line emission to the nucleosynthesis. For this purpose we have for many of 
the lines compared the llEl model by Nomoto et al. and the lOH model by Woosley. First we 
note that the most apparent discrepancy is in [CI] A A 9825, 9850, which is similar in both 
models. This over-production, together with the observed presence of [N I] A A 10398, 10408 lines, 
indicates a less extensive zone of He - C, replaced by a He - N zone of CNO-burning products. 
This illustrates that although different models agree, there may be systematic effects due to 
uncertainties in the input physics in the models. 

The main differences between the two models are found in the Si - S, and oxygen-rich zones. 
In the llEl model, neon and magnesium are, next to oxygen, the most abundant elements, while 
silicon and sulphur replaces these in the lOH model. The effects of this are most prominent in the 
light curves of the [0 1] A A 6300,6364, [Ne II] A 12.814 /im, Mg I] A 4571, [Si I] A 1.6454 ^um, 
and [S I] A 1.0820 /im hnes. 

Starting with oxygen, the thermal part of the light curve is well reproduced in both models, 
although the llEl model gives a somewhat better fit. When it comes to the non-thermal part, 
both models under-produce the light curve by a large factor, although lOH somewhat less than 
llEl. The problem with the [O I] line is discussed further below. 

[Ne II] A 12.814 fim clearly favors the lOH model, because of the smaller neon mass in the 
oxygen zone, while Mg I] A 4571 agrees on the average better in the lOH model. This applies 
especially to the lOH model where we decreased the photoionization rate to zero (section [5.31 ). In 
the llEl model the light curve has the right form, but is at all epochs a factor of 2 - 3 to strong. 

Because of the dominance of the zones outside of the Si - S zone for the [Ca II] emission, 
these lines are not as sensitive to the differences in the Si - S zone between the two models, as one 
could otherwise have hoped. Because of the larger Ca II contribution from the O - Si - S zone in 
the lOH model earlier than 400 days, this model gives a higher flux and better agreement with the 
observations compared to the llEl model. 

While the contribution from the Si - S zone to the [Si I] - [Fe II] A 1.64 /xm blend varies 
considerably, the [Si I] emission is dominated by the oxygen zone and does not differ much. The 
[S I] A 1.0820 /im line is in lOH over-produced by a factor of ~ 2 at 200 - 300 days, and favors 
therefore the llEl model. 
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Summarizing this evaluation, there is no clear winner. While the [S I] lines favor llEl, [Ne 
II], Mg I] and Ca II lines agree better with lOH. The early evolution of the important [O I] lines 
may favor the llEl model, but because of the serious discrepancy at late time, this is not obvious. 
Most discrepancies are at the factor of two level, and in view of other uncertainties it is not clear 
how significant the differences are. One may also take the opposite view and conclude that within 
this factor there is good agreement with most abundances in the models. 

The low level of the non-thermal plateau of the [0 1] A A 6300, 6364 lines poses a problem for 
our 'standard' models. We have already discussed various ingredients in a solution to this problem. 
From our analytical discussion it is obvious that an appreciable contribution to the excitation 
from the triplets is needed. Although the energy to the triplets is sufficient, and branching to the 
singlets is probable, photoabsorption of the A 1302 line by Mg I and Si I decreases this mode by 
a large factor. As Figure ^ shows, only a decrease in the photoabsorption to near zero has an 
appreciable effect on the non-thermal level. Without a strong decrease of both the magnesium 
and silicon abundance this is difficult to achieve, unless we have under-estimated the ionization 
of Mg I and Si I by a large factor. In the explosion models either magnesium or silicon, or both, 
tend to have a high abundance in the inner oxygen region. It is somewhat ironic that the most 
simplified model, with only oxygen, is the most successful model for the non-thermal phase. We 
have, however, also seen that the maximum triplet contribution from recombination in most cases 
is accompanied by too strong O I recombination lines, and we have argued that charge transfer 
with Si I may be a solution to this. 

Another possible solution to the [O I] problem would be that the electron fraction is much 
lower than we have calculated, because the excitation efficiency increases with decreasing electron 
fraction. We have already experimented with a higher density (i.e. lower ffiling factor), and this 
did not have much effect, because Xg oc n~^/^. Charge transfer can probably only change by 
factors of a few, because the recombination rates of Mg II, Si II, and Ne II are similar to that of 
O II. The UV field does not influence the O I emission appreciably. 

A lower abundance of silicon in the O - Si - S region would help improving both the [Si I] 
A 1.6454 /im line, and also increase the non-thermal plateau of the [0 1] A A 6300,6364 lines. In 
fact, a general decrease of the convective mixing of elements, like silicon in the O - Si - S region 
and carbon in the He - C region, in connection with the shell burning stages, would probably 
improve the agreement with the observations considerably. It is interesting that more recent 
explosion models by Woosley & Weaver (1995) give lower silicon abundances. Their 15 Mq model 
has X{Si) ~ 1.2 x 10"^ and the 25 M© model X(Si) ~ 3.8 x 10'^, both lower than in the llEl 
model. These models also show that the extent of the He - C zone can vary substantially. In 
the 15 Mq model the He - C zone accounts for less than half of the helium zone, with the rest 
is He - N. In the 25 Mq model nearly all of the helium zone is in the form of He - C. The neon 
abundance is in these models higher than even in the llEl model, possibly creating problems. 

We also remark that the exact form of the continuum destruction probability is important. 
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and unfortunately uncertain. This is illustrated by the difference between equations (32), (33) and 
(34) in Paper I, where /3c varies by factors of ten or more, depending on the value of kc/ki^- A 
further study of this for the conditions relevant for the ejecta would be of great interest. 

Because most published observations were obtained during the first four years, we have in 
this paper concentrated on the evolution up to ~ 1200 days. Our models are, however, calculated 
up to 2000 days, and we therefore briefly compare these with the HST observations by Wang et 
al. (1996) at 1862 days. At this epoch blending in the UV is less severe than at earlier epochs. 
Meaningful estimates of the fluxes of e.g. Mg I A 2852 and Mg II A 2800, as well as the Balmer 
continuum, can therefore be obtained. From a fit to the Balmer continuum Wang et al. find 
a temperature at 1862 days of 500 it 100 K. The temperature we find in the hydrogen-rich gas 
decreases from ~ 425 K in the core to ~ 390 K in the inner envelope and ~ 220 K at the highest 
velocities (Figures 2 and 7 in Paper I). Given the observational uncertainties, especially in the 
level of the 'continuum', we find that this is in reasonable agreement with the observations. 

In Table |^ we give the observed fluxes at 1862 days from Wang et al. compared to the 
reddening adjusted fluxes from our lOH model at the same epoch. We find the Balmer continuum 
flux to be a factor of two larger than observed. The latter may, however, be a lower limit because 
of the high level of the background flux assumed. A factor of two higher fiux could well be 
consistent with the observations. Our Ha fiux is very close to the observed value, implying that 
our electron density is close to that derived from the observations. Wang et al. estimate from Hct 
an electron density of ~ 6 x 10^ cm~^, while we find 4.4 x 10^ — 1.6 x 10^ cm~^ in the different 
regions. The [0 1] A A 6300, 6364 flux is in the model a factor of two lower than the observed. 
Although the discrepancy is smaller than at earlier epochs, there is still a substantial uncertainty 
in the model flux because of the effects discussed earlier in this section. The Mg I A 2852 and 
Mg II A A 2796, 2803 lines are not resolved in the observations. The total flux of the blend is a 
factor of ~ 2.6 higher than the calculated value. In the future we hope to make a more detailed 
comparison with these observations, as well as those even later in Chugai et al. (1997). 

7. CONCLUSIONS 

For convenience we here summarize our main results. 

1. From the evolution of especially the [O I] and [Fe II] lines there is overwhelming evidence 
that the IR-catastrophe has occurred in the metal-rich regions at ~ 700 days. 

2. The nearly constant ratio of the [O I] A 6300, 6364 luminosity and the bolometric luminosity 
at t ^ 800 days shows that non-thermal excitation dominates this part of the light curve, 
which implies a very low temperature. The same applies to e.g. Ha, He I A 2.058 ;um, and 
Mg I] A 4571, although here recombination following non-thermal ionization dominates. 



- 35 - 



For reliable mass estimates, constraints from the line profiles have to be taken into account. 
Even faint line wings at high velocity may correspond to comparatively large masses. The 
most reliable masses are therefore those within ~ 3000 km s~^. 

We find that the density in the hydrogen envelope is considerably flatter than in published 
models, decreasing like p oc V~'^, out to ~ 5000 km s~^. The total hydrogen zone mass is 
~ 7.7 Mq, of which ~ 3.9 is hydrogen. The mass of the hydrogen-rich gas within 2000 
©• 



km s"^ is ~ 2.2 M, 



5. The helium mass is ~ 5.8 Mq, of which ~ 1.9 is in the helium zone and ~ 3.9 Mq in 
the hydrogen zone. Although the He I A 2.058 /im line profile is uncertain, models and 
observations indicate that the helium emitting mass may extend to ~ 4000 km s~^. In the 
center it penetrates to <^ 500 km s~^. Before ~ 500 days the light curve, and therefore the 
mass estimate, is sensitive to assumptions about complete versus partial redistribution. 

6. Our best determinations of the other masses are M(N) ~ 3.4 x 10^^ Mq, M(Ne) 

6 X 10-2 Mq, M(Mg) « 2.2 x lO^^ Mq, M{^^m + ^Om) « 6 x lO'^ Mq. The total masses 
of carbon, sulphur and calcium are not well determined, because of the low contributions 
from zones where these elements are most abundant. 

7. The oxygen distribution extends to <^ 400 km s^"*^, and is consistent with ~ 1.9 Mq of 
oxygen enriched gas. Because of uncertainties connected with the non-thermal part of the 



light curve there could be a considerable uncertainty in this number, although Figure 11 
gives the likely range. The [O I] line therefore remains an unsolved problem and a better 
understanding of this line is clearly needed. 

8. The extent of convective mixing of carbon into the helium-rich region is likely to be 
more limited than stellar evolution calculations show. Instead, a more massive zone of 
CNO-burning products is indicated by the observations. 

9. Most lines, including the O I, Ca II, and Fe II lines, have contributions from more than 
one component. Even trace amounts in zones with a large gamma-ray deposition may give 
substantial contributions to the total line fluxes. 

10. The [Fe II] lines have a large contribution from primordial iron, and at late epochs this 
component dominates, because of the IR-catastrophe in the iron core. The extended wings 
of the [Fe II] lines are probably produced by primordial, rather than synthesized, iron. 

11. Both charge transfer and photoionization effects are important for the strengths of lines from 
trace ions. The observations indicate strong suppression of the UV flux by dust absorption 
and resonance scattering. 



Finally, after this paper was submitted we received a preprint by de Kool, Li, & McCray 
(1997) where a similar analysis and calculation is performed. In spite of detailed differences, it is 
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comforting that the results of these two independent calculations with regard to both the physical 
conditions and the observed emission are very similar. 
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Table 1. Modeled and observed HST-fluxes at 



1862 days. 


Line 


ModeP 


Observation*^ 


Ha 


18.8 


16.6 


Balmer cont. 


7.2 


3.1 


[0 I] 6300, 6364 


1.2 


2.9 


Mg I] 4571 


4.6 


2.2 


Mg I 2852 


0.58 


1.7^ 


Mg II 2796, 2803 


0.07 





^Observed fluxes are from Wang et al. (1996), 
while the modeled are from the lOH model, 
adjusted for reddening. All fluxes are in units of 
10^^^ erg cm~^. 

^'Total flux for Mg I A 2852 and Mg II 
A A 2795, 2803. 



-38- 



REFERENCES 

Aitken, D. K., Smith, C. H., James, S. D., Roche, P. F., Hyland, A. R., & McGregor, R J. 1988, 
MNRAS, 235, 19p 

Arnett, W. D. 1996, Supernovae and Nucleosynthesis, Princeton University Press, Princeton 
Basko, M. 1994, ApJ, 425, 264 

Bautista, M. A., DePoy, D. L., Pradhan, A. K., Ehas, J. H., Gregory, B., Philhps, M. M., k 
Suntzeff, N. B. 1995, AJ, 109, 729 

Bouchet, P., & Danziger, I. J. 1993, A&A, 273, 451 

Bouchet, P., Danziger, I. J & Lucy, L. B. 1991, AJ, 102, 1135 

Bouchet, P., Danziger, I. J., Gouiffes, C., DeUa Valle, M., k Moneti, A. 1996, in lAU. Coll. 

No. 145, Supernovae and Supernova Remnants, eds. R. McCray &; Z. Wang, Cambridge 
University Press, 201 

Chugai, N. N. 1987, Astrophysics, 26, 53 

Chugai, N. N., Chevalier, R. A., Kirshner, R. P., & Challis, P. M. 1997, ApJ, 483, 925 
Chung S., Lin, C. C, & Lee, E. T. P. 1991, Phys. Rev. A, 43, 3433 

Colgan, S. W. J., Haas, M. R., Erickson, E. F., Lord, S. D., & HoUenbach, D. J. 1994, ApJ, 427, 
874 

Danziger, I. J., Bouchet, P., Gouiffes, C, & Lucy, L. B. 1991, in Proc. ESO/EIPC Supernova 
Workshop, SN 1987A and other Supernovae, eds. LJ. Danziger & K. Kjar, (Garching: 
ESQ), 217 

Dwek, E., Moseley, S.H., Glaccum, W., Graham, J.R., Loewenstein, R.F., Silverberg, R.F., &; 
Smith, R. K. 1992, ApJ, 389, L21 

Erickson, E. F., Haas, M. R., Colgan, S. W. J., Lord, S. D., Burton, M. G., Wolf, J., Hollenbach, 
D. J., & Werner, M. 1988, ApJ, 330, L39 

Ferland, G. J. & Persson, S. E. 1989, ApJ, 347, 656 

Fitzpatrick, E. L. 1985, ApJ, 299, 219 

Fransson, C, 1984, A&A, 132, 115 

Fransson, C, 1994, in Supernovae (Les Houches, Session LIV 1990), eds. J. Audouze, S. Bludman, 
R. Mochkovitch, & J. Zinn- Justin (New York: Elsevier), 677 

Fransson, C, & Chevalier, R. A. 1989, ApJ, 343, 323 

Fransson, C, Houck, J. Sz Kozma, C. 1996, in lAU. Coll. No. 145, Supernovae and Supernova 
Remnants, eds. R. McCray &; Z. Wang, Cambridge University Press, 211 

Fransson, C, Cassatella, A., Gilmozzi, R., Kirshner, R. P., Panagia, N., Sonneborn, G., & 
Wamsteker, W. 1989, ApJ, 336, 429 



-39- 



Pryxell, B., Miiller, E., & Arnett, D. 1991, ApJ, 367, 619 

Haas, M. R., Colgan, S. W. J., Erickson, E. F., Lord, S. D., Burton, M. G., & Hollenbach, D. J. 
1990, ApJ, 360, 257 

Hashimoto, M., Nomoto, K., & Shigeyama, T. 1989, A&A, 210, L5 

Herant, M., & Benz, W. 1992, ApJ, 387, 294 

Houck, J. C, & Fransson, C, 1996, ApJ, 456, 811 

Hummer, D. G. 1968, MNRAS, 138, 73 

Hummer, D. G., & Rybicki, G. B. 1985, ApJ, 293, 258 

Jennings, D. E., Boyle, R. J., Weidemann, G. R., & Moseley, S. H. 1993, ApJ, 408, 277 

Julienne, RS., Davies, J. &; Oran, E. 1974, J. Geophys. Res. 79, 2540 

Kozasa, T., Hasegawa, H., & Nomoto, K. 1991, A&A, 249, 474 

Kozma, C, & Fransson, C. 1992, ApJ, 390, 602 (KF92) 

Kozma, C, & Fransson, C. 1997, ApJ, previous paper (Paper I) 

Li, H., & McCray, R. 1992, ApJ, 387, 309 

Li, H., & McCray, R. 1993, ApJ, 405, 730 

Li, H., & McCray, R. 1995, ApJ, 441, 821 

Li, H., & McCray, R. 1996, ApJ, 456, 370 

Li, H., McCray, R., & Sunyaev, R. A. 1993, ApJ, 419, 824 

Liu, W. & Dalgarno, A. 1995, ApJ, 454, 472 

Lucy, L. B., Danziger, I. J., Gouiffes, C, & Bouchet, P. 1989, in lAU Colloquium No. 120, 

Structure and Dynamics of the Interstellar Medium, eds. G. Tenorio-Tagle, M. Moles, and 
J. Melnick, Lecture Notes in Physics (Springer Verlag), 164. 

Lucy, L.B., Danziger, I.J., Gouiffes, C, &: Bouchet, P. 1991, in Supernovae, Proc. of the Tenth 
Santa Cruz Summer Workshop in Astronomy and Astrophysics, ed. S.E. Woosley (Springer 
Verlag), 82. 

Martin, P. G. 1988, ApJS, 66, 125 

McGregor, P. J. 1988, Proc. Astron. Soc. Austraha, 7, 45 

Meikle, W. P. S., Allen, D. A., Spyromilio, J., & Varani, G. -F. 1989, MNRAS, 238, 193 

Meikle, W. P. S., Spyromilio, J., Allen, D. A., Varani, G. -F., & Gumming, R. J. 1993, MNRAS, 
261, 535 

Menzies, J.W. 1991, in Proc. ESO/EIPC Supernova Workshop, SN 1987A and other Supernovae, 
eds. LJ. Danziger & K. Kjar, (Garching: ESQ), 209 

Moseley, S. H., Dwek, E., Glaccum, W., Graham, J. R., Loewenstein, R. F., & Silverberg, R. F. 
1989, ApJ, 347, 1119 



-40- 



Nomoto, K., & Hashimoto, M. 1988, Phys. Rept. 163, 13 
Oliva, E. 1992, priv. comm. in Li, McCray, & Sunyaev (1993) 
Oliva, E. 1993, A&A, 276, 415 

Oliva, E., Moorwood, A. F. M., & Danziger, I. J. 1989, in 22nd ESLAB symp.. Infrared 
Spectroscopy in Astronomy, ed. B. H. Kaldeich, ESA SP-290, 375 

Phillips, M. M., Hamuy, M., Heathcote, S. R., SuntzefF, N. B., & Kirhakos, S. 1990, AJ, 99, 1133 

Phillips, M. M., & Wilhams, R. E. 1991 in Supernovae, ed. S. E. Woosley (New York: Springer), 
36 

Pun, C. S. J., Kirshner, R. P., Sonneborn, G., Challis, P., Nassiopoulos, G., Arquilla, R., 
Crenshaw, D. M., Shrader, C, Teays, T., Cassatella, A., Gilmozzi, R., Talavera, A., 
Wamsteker, W., Fransson, C., & Panagia, N. 1995, ApJS, 99, 223 

Roche, P. F., Aitken, D. K., k Smith, C. H. 1993, MNRAS, 261, 522 

Savage, B. D., & Mathis, J. 1979, ARA&A, 17, 73 

Shigeyama, T., & Nomoto, K. 1990, ApJ, 360, 242 

Shigeyama, T., Nomoto, K., & Hashimoto, M. 1988, A&A, 196, 141 

Shull, J. M., & Van Steenberg, M. E. 1982, ApJS, 48, 95 

Sonneborn, G., Fransson, C., Lundqvist, P., Cassatella, A., Gilmozzi, R., Kirshner, R. P., Panagia, 
N., & Wamsteker, W. 1997, ApJ, 477, 848 

Spyromilio, J., & Graham, J. R. 1992, MNRAS, 255, 671 

Spyromilio, J., & Pinto, P. A. 1991, in Proc. ESO/EIPC Supernova Workshop, SN 1987A and 
other Supernovae, eds. I.J. Danziger &; K. Kjar, (Garching: ESO), 423 

Spyromilio, J., Meikle, W. P. S., & Allen, D. A. 1990, MNRAS, 242, 669 

Spyromilio, J., Stathakis, R. A., Cannon, R. D., Waterman, L. & Couch, W.J. 1991, MNRAS, 
248, 465 

Storey, P. J., & Hummer, D. G. 1995, MNRAS, 272, 41 

SuntzefF, N. B., Phillips, M. M., Depoy, D. L., Elias J. H., & Walker, A. R. 1991, AJ, 102, 1118 
Swartz, D. A. 1994, ApJ, 428, 267 

Thielemann, F.-K., Nomoto, K., & Hashimoto, M. 1996, ApJ, 460, 408 

Varani, G.-F., Meikle, W. P. S., Spyromilio, J., & Allen, D.A. 1990, MNRAS, 245, 570 

Wang, L., Wheeler, J. C, Kirshner, R. P., Challis, P. M., Filippenko, A. V., Fransson, C, 
Panagia, N., Phillips, M. M., & SuntzefF, N. 1996, ApJ, 466, 998 

Whitelock, P. A. et al. 1989, MNRAS, 240, 7P 

Wooden, D. H., Rank, D. M., Bregman J. D., Witteborn, F. C, Tielens, A. G. G. M., Cohen, M., 
Pinto, P. A., & Axelrod, T. S. 1993, ApJS, 88, 477 



-41 - 



Woods, D. T., ShuU, J. M., & Sarazin, C. L. 1981, ApJ, 249, 399 
Woosley, S. E. 1988, ApJ, 330, 218 

Woosley, S. E., Pinto, P. A., & Hartmann, D. 1989, ApJ, 346, 395 
Woosley, S. E., & Weaver, T. A. 1986, ARA&A, 24, 205 
Woosley, S. E., & Weaver, T. A. 1995, ApJS, 101, 181 
Xu, Y., McCray, R., Oliva, E., & Randich, S. 1992, ApJ, 386, 181 



This preprint was prepared with the AAS lAT^^X macros v4.0. 



-42- 




Velocity (km s ^) 

Fig. 1. — The Ha line profile at ~ 800 days. The observed line profile from Phillips et al. 
(1990) (solid line), compared to two of our models for the hydrogen distribution. The bump 
at ~ 933 km is due to an unresolved, narrow [N II] A 6583 line from the ring. Note the too 
weak line wings and the too fiat peak in the 14E1 model. 
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Fig. 2. — The time evolution for different hydrogen hnes (sohd hnes), together with observations 
(open dots). Data from Danziger et al. (1991), Bouchet &; Danziger (1993), Wooden et al. (1993), 
Meikle et al. (1989, 1993), and Bautista et al. (1995). The contributions from the core region are 
shown as dotted Hnes and from the envelope as dashed lines. 
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Fig. 3. Same as Figure 2 but for the Balmer and Paschcn continuum, as well as for the H 9— >7, 
and the H 7^6 hnes. Observations are from Roche, Aitken, k, Smith (1993), and Wooden et al. 
(1993). 
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Fig. 4. — Time evolution of the He I A A 1.0830 /xm and 2.058 /xm lines. Calculations for both the 

lOH and llEl models are shown together with observations. In the two upper figures the sum of 

the He I A 1.0830 ^m, [S I] A 1.0820 //m and Pa7 are shown, together with the contributions from 

the different composition zones. In the middle figures the contributions from all zones due to He I 
\ 1 noon ..™ TQ Tl \ i neon i^u^ f;™,„«.^ ^u^,„ j-U« +^+.,1 
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Fig. 5.— The He I A 2.058 /xm line at 574 and 695 days from Meikle et al. (1993), together 
with our calculated profile at 650 days (solid line), including dust. The helium distribution used is 
discussed in the text. 
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Fig. 6. — The effect of different continuum destruction probabilities for He I A 2.058 /xm in tfie 
lOH model. The solid lines show the total luminosity, including all zones, while the dotted lines 
show the contribution from the He - C zone alone. The upper curves, labeled HR85, give the light 
curves for the destruction probability given by Hummer & Rybicki (1985) for a Doppler profile, 
assuming complete redistribution, while the lower curves, labeled H68, gives the same for a Voigt 
profile and complete redistribution from Hummer (1968). 
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Fig. 7. — Light curve of [CI] A A 9824,9850 for the lOH and llEl models. Observations are from 
SpyromiHo et al. 1991. 
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Fig. 8.— Light curves of [CI] A A 9824, 9850, [N I] A A 10398, 10408 and [N II] A A 6548, 6583 
for a model where the He - C zone in the lOH model has been replaced by a He - N zone. 




Fig. 9.— 



Grotrian diagram of the levels and transitions included in our O I model atom. 
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Fig. 10. — Light curve of [0 1] A A 6300, 6364. The upper figure shows the fluxes from the lOH 
model and the lower from the llEl model. Observations are from Danziger et al. (1991). 
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Fig. 11. — [0 1] A A 6300, 6364 light curves for three different masses of the oxygen zones 
in the llEl model. The solid line shows the standard model with M(0) = 1.9 Mq. Note that 
while the thermal part of the light curve is probably bracketted by this range, the non-thermal is 
under-produced even for the highest mass. 
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Fig. 12. — [0 1] A A 6300, 6364 light curves for the llEl model, showing the importance of 
the photoelectric absorption by Si I and Mg I. The solid line is the standard model for the O 
- Ne - Mg region with X(Si) = 1.64 x 10~^, while the dotted line is a model with the silicon 
abundance decreased to 1.3 x 10~^. Note the increase in the non-thermal plateau at times later 
than ~ 900 days, and the bump between 600 - 900 days. The dashed line shows the flux for a model 
where the photoelectric absorption has been switched off for both Si I and Mg I. This increases the 
non-thermal flux by a factor three, compared to the standard model. 



-54- 



[0 I] 6300, 6364 

"1 \ \ \ \ \ \ \ \ \ \ \ \ \ \ [ 




J \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ L 

-6000 -4000 -2000 2000 4000 6000 

Velocity (km s ^) 



Fig. 13.^ The [0 1] A A 6300, 6364 line profile for the llEl raodel at 800 days (dotted line), 
together with observations from Phillips et al. (1990) (solid line). The more peaked observed profile 
shows that mixing of oxygen reaches further in, to less than 400 km s~^, than in the model. 
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Fig. 14. — Light curve of [Ne II] A 12.814 for the lOH and the llEl model, respectively. 
Observations are from Aitken et al. (1988), Roche et al. (1993), and Colgan et al. (1994). 
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Fig. 15. — Light curve of A Mg I] 4571 for the lOH and the llEl model. Observations are from 
Danziger et al. (1991). 
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Fig. 16.— Light curve of Mg II A A 2795, 2802 for the lOH and the llEl model. 
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Fig. 17. — The upper two figures show the contributions from Fe II and Si I to the 1.64 /im feature 
for the lOH and llEl model, respectively. The lower two figures show the contribution to the [Si 
I] A 1.6454 fim fine from the different composition zones for the two models. Observations from 
Meikle et al. (1989), (1993), and Bautista et al. (1995). 
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Fig. 18.— The upper figures show the fight curve of [Ca II] A A 7291, 7324 for the lOH and llEl 

models, respectively, with contributions from the individual components. In the lower figures the 
light curve of the [Ca II] A A 7291, 7324 lines, as well as for other strong lines in the wavelength 
range 7200 - 7400 A are shown. The curve labeled Fe II includes several weak lines. Observations 
from Spyromilio et al. (1991), from Phillips &; Williams (1991) and from Suntzeff et al. (1991). 
Note the strong deficiency later than ~ 800 days, probably indicative of photoexcitation. 
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Fig. 19. — The upper figures show the Ca II IR-triplet with contributions from the individual 
regions for lOH and llEl, respectively. In the lower figures the light curve of other lines in the 
wavelength range 8400 - 8750 A are shown. The curve labeled H I contains the Paschen lines in 
this wavelength region. Observations from Spyromilio et al. (1991) and from Suntzeff et al. (1991). 
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Fig. 20. — Light curves for some of the most important Fe II lines. The calculations based on the 
lOH and llEl models do not differ significantly from each other. The contributions from primordial 
and processed iron are shown separately. Because of the IR-catastrophe the [Fe II] flux from the 
iron core drops to a very low level at t ;^ 600 days, except for the A 25.99 /xm line. References for 
the observations are given in the text. 
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Fig. 21.— Light curves of [Co II] A 10.52 /urn, [Ni II] A 6.634 /xm, and [Ni II] A 10.68 /xm. 
The calculations based on the lOH and llEl models do not differ significantly from each other. 
Observations are from Aitken et al. (1988), Jennings et al. (1993), Roche et al. (1993), Wooden et 
al. (1993), Colgan et al. (1994). 
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Fig. 22.— Light curves of [Fe I] A 1.44 //m, [Ni I] A 3.119 fim, and [Ni I] A 7.505 fim. Note that 

in modeling these lines no photoionization of these elements is included. The calculations based on 
the lOH and llEl models do not differ significantly from each other. Observations for [Fe I] A 
1.44 /im are taken from Oliva (1992) and Oliva, Moorwood, & Danziger (1989). Observations for 
the [Ni I] lines are from Meikle et al. (1989, 1993), and Wooden et al. (1993). 



